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Abstract. The classical one-phase Stefan problem describes the temperature distribution in a 
homogeneous medium undergoing a phase transition, such as ice melting to water. This is accom- 
plished by solving the heat equation on a time-dependent domain whose boundary is transported by 
the normal derivative of the temperature along the evolving and a priori unknown free-boundary. 
We establish a global-in-time stability result for nearly spherical geometries and small tempera- 
tures, using a novel hybrid methodology, which combines energy estimates, decay estimates, and 
Hopf-typc inequalities. 
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1. Introduction 

1.1. The problem formulation. We consider the problem of global existence and asymptotic 
stability of classical solutions to the classical Stefan problem describing the evolving free-boundary 
between the liquid and solid phases. The temperature of the liquid p(t,x) and the a priori unknown 
moving phase boundary T(t) must satisfy the following system of equations: 

p t -Ap = in fi(i); (1.1a) 

V(T(t)) = -d n p on r(f); (1.1b) 

p = on T(t); (1.1c) 

p(o,-)=po,fi(o)=no. (l.id) 

For each instant of time ts[0,T], f2(i) is a time-dependent open subset of K d with d>2, and 
r(t) :=dfl(t) denotes the moving, time-dependent free-boundary. 

The heat equation (ll.lal) models thermal diffusion in the bulk f2(t) with thermal diffusivity set 
to 1. The boundary transport equation (jl.lbl) states that each point on the moving boundary is 
transported with normal velocity equal to —d n p = — Vp-n, the normal derivative of p on T(t). Here 
n denotes the outward pointing unit normal to T(t), and V(T(t)) denotes the speed or the normal 
velocity of the hypersurface T(t). The homogeneous Dirichlct boundary condition (|l.lc|) is termed 
the classical Stefan condition and problem (|l.ljl is called the classical Stefan problem. It implies that 
the freezing of the liquid occurs at a constant temperature p = 0. Finally, we must specify the initial 
temperature distribution po:flo— >R, as well as the initial geometry Sin. Because the liquid phase 
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is characterized by the set : p(x,t)>0}, we shall consider initial data po > in no- 

Problem belongs to the category of free boundary problems which are of parabolic-hyperbolic 
type. Thanks to (|l.lal) . the parabolic Hopf lemma implies that <9„p(t)<0 on T(t) for t>0, so we 
impose the non- degeneracy condition or so-called Taylor sign conditioi\j 

~d nPo >X>0 on T(0) (1.2) 

on our initial temperature distribution. Under the above assumptions, we proved in Hadzic & 
Shkoller [29] that ([1-ljl is indeed well-posed. 

1.2. The reference domain fl and the dimension. For our reference domain, we choose the 
unit ball in R 2 given by 

fi = B(0,l):={a:eR 2 : M<1}, 

with boundary r = S 1 :={i6l 2 : |x| = 1}. We shall consider initial domains fio whose boundary Tq 
is a graph over the reference boundary T. In order to simplify our presentation, we consider evolving 
domains Q(t) in R 2 , but as we shall explain in Section^ our methodology works equally well in any 
dimension d> 2. 

Our choice of the reference domain fl follows from two considerations. First, we need employ 
only one global coordinate system near the boundary T, rather than a collection of local coordinate 
charts that a more general domain would necessitate, and the use of one coordinate system greatly 
simplifies the presentation of our energy identities, that provide very natural estimates for the second- 
fundamental form of the evolving free-boundary T(t). Second, we shall need quantitative Hopf- type 
inequalities in order to bound the term defined in f| 1 . 2|) from below, and such estimates are available 
in a particularly satisfying form in the case of the nearly spherical domains, thanks to the explicit 
construction of comparison functions in Oddson [40] . 

1.3. Notation. For any s > and given functions / : 0— > R, ^:T->Kwe set 

\\f\\s--=\\f\\H°(n) and \ip\ s := ||^|| ff .(r)- 

H s (n)' shall denote the dual space of while on the boundary, H s (r)' = H- S (T). If i = 1,2 

then f,i'-=d x if is the partial derivative of / with respect to x l . Similarly, f,ij'=d x id x jf, etc. For 
time-differentiation, f t :=dtf '■ Furthermore, for a function f(t,x), we shall often write f(t) for f(t,-), 
and /(0) to mean f(0,x). We use B:=t- V to denote the tangential derivative, so that 

8f:=d e f, B k f: ;)!;/. 

where 0£ [0,2ir) denotes the angular component in polar coordinates. The Greek letter a will often 
be reserved for multi-indices a = («i,a2), with d a :=d x *d x ^ and \a\ =ai +«2- The identity map on 
n is denoted by e{x) =x, while the identity matrix is denoted by Id. We use C to denote a universal 
(or generic) constant that may change from inequality to inequality. We write X < Y to denote 
X<CY. We use the notation P(-) to denote generic real polynomial function of its argument(s) 
with positive coefficients. The Einstein summation convention is employed, indicating summation 
over repeated indices. The L 2 -inner product on fi is denoted by (-,-)i2. 

1.4. Fixing the domain. We transform the Stefan problem (|1.1[) . set on the moving domain £l(t), to 
an equivalent problem on the fixed domain £1. For many problems in fluid dynamics, the Lagrangian 
flow map of the fluid velocity provides a natural family of diffcomorphisms which can be used to fix 
the domain, but for the classical Stefan problem, we use instead (in the parlance of fluid dynamics) 
the so-called Arbitrary Lagrangian-Eulerian (ALE) family of diffcomorphisms; these ALE maps 
interpolate between the Lagrangian and Eulerian representations of the equations. For this problem, 
we choose a simple type of ALE map, consisting of harmonic coordinates, also known as the harmonic 
gauge. 



^This type of stability condition dates back to the early work of Lord Rayleigh 1441 and Taylor 1461 in fluid 
mechanics, and appears as a necessary well-poscdness condition on the initial data in many free-boundary problems 
wherein the effects of surface tension are ignored; examples include the Helc-Shaw cell, the water waves equations 
|48| , and the full Euler equations in both incompressible |15| and compressible form |18l I17| . 
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1.4.1. The diffeormorphism 4*(i). We represent our moving domain fl(t) as the image of a time- 
dependent family of diffeomorphisms ^(t) :f2— >Q,(t). In order to define these diffcomorphisms, we 
let h(t, •) : r — denote the signed height function whose graph (over T) is the set T(t). For £ € T = S 1 , 
we define the map 

*(t,0 = (i+M*»OK=-R(*»OS 

which is a diffcomorphism of T onto T(t) as long as h(t) remains a graph. The outward-pointing 
unit normal vector n(t,-) to the moving surface T(t) is defined by 

(nom,o=mi/\mi\- 

We shall henceforth drop the explicit composition with the diffeomorphism and simply write 

n(t,o=mi/\(mi\ 

for the unit normal to the moving boundary at the point ^(t,^) £T(t). 

Introducing the unit normal and tangent vectors to the reference surface T as 

N:=£, T-.= £ g or equivalently N(6) = (cos0,sin0), t(6) = (-sine*, cosfl) , (1.3) 

we write the unit normal to T(t) as 

. ^ RN — hgT , ,, 

n (t£) = — (1.4) 

The evolution of h(t) is then given by 

h t = vN(9)-^v-T(6). (1.5) 
H 

Assuming that the signed height function h(t,-) is sufficiently regular and remains a graph, we 
can define a diffeomorphism ^:f2— >Q(t) as the elliptic extension of the boundary diffeomorphism 
£ i-> (l + h(£,t))£, by solving the following Dirichlet problem 

A* = in Q, 

*(t,0 = -R(U)£ Cer. (1.6) 

Since the identity map e : fi — > is harmonic in ft and "J — e = h£ on T, standard elliptic regularity 
theory for solutions to (|1.6p shows that 

\\*-e\\ H -w<C\\h\\ H .- . Hr) , s>0.5, (1.7) 

so that for h(t) sufficiently small and s large enough, the Sobolev embedding theorem shows that 
V\P(t) is close to the identity matrix Id, and by the inverse function theorem, each ^(i) is a diffco- 
morphism. 

1.4.2. The temperature and velocity variables on the fixed domain Q. First we introduce the velocity 
variable u = — Vp in Q(t). Next, we introduce the following new variables set on the fixed domain 

n-. 

q = poty (temperature), 

v = uoty (velocity), 

w = ^ t (extension of boundary velocity vector), 

A=[D n< l>\~ 1 (inverse of the deformation tensor), 

J = detD^ ( J acobian determinant), 

a=JA (cofactor matrix of the deformation tensor). 

The relation u = — Vp is then written as v l + A^q,k — for i = 1,2. By the chain rule, 

q t —p t o * + (\7po *) ■ * t =p t o - v ■ w , 

and 



ApoV = Ayq:=Ai(A l ?q, k ), 



j ■ 
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Letting n = J^ 1 (R 2 + R 2 )^n, we see that 

h i {t,x)=A k i {t,x)N k {x), (1.8) 

and equation (|1.5p can be written as h t = v-n/Rj, where Rj = RJ^ 1 . Note that Rj = RJ^ 1 = 
(l + h),!^ 1 is very close to 1. 

1.4.3. The classical Stefan problem set on the fixed domain CI. The classical Stefan problem on the 
fixed domain Cl is written as 

qt-Ai(Afq, k ) tj = -v* t in (0,T]xfi, (1.9a) 

v i +A k q, k =0 in [0,T]xO, (1.9b) 

q = on [0,T]xI\ (1.9c) 

h t = v-N-(l + h)- 1 h e v-T on (0,T]xl, (1.9d) 

A* = on [0,T]xQ, (1.9c) 

* = (l + /i)iV on [0,T]xr, (1.9f) 

q = q o >0 on {t = 0}xft, (1.9g) 

ft = ft on {f = 0}xr, (1.9h) 

where the initial boundary OCIq is given as a graph over Cl with the initial height function ho, i.e. 
afi = {ieR 2 , z = (l + fto(0)£) C^S 1 }. Note that $ = <5(0) : Cl^ ft is a near identity transforma- 
tion, mapping the reference domain Cl onto the initial domain f^o- The initial temperature function 
go equals po ^- Problem (|1.9j) is a reformulation of the problem (|1.1[) . 

Henceforth, without loss of generality, wc shall assume that the initial domain CIq is the unit ball 
-Bi(O) or in other words ho = 0. In this case, we set $ = e, where e : £1 — ^ S7 is the identity map, and 
^{t)\t=o = e - I n Section [5l we will explain the minor modification required when /i ^0, as well as 
the case that the dimension d~3. 

Observe that the boundary condition (|1 .9d|) implies that 

ty t -n(t)=v-n(t) on [0,T]xT so that *(i)(r) =T(t) . (1.10) 

1.4.4. The energy and dissipation functions. Near T = dCl, it is convenient to use tangential deriva- 
tives d:=de with 8 denoting the polar angle, while near the origin, Cartesian partial derivatives d Xi 
are natural. For this reason, we introduce a non-negative C°° cut-off function fi:Cl—>M.+ with the 
property 

fx{x) = if \x\<l/2; n{x) = \ if 3/4< \x\ < 1. 

Definition 1.1 (Higher-order norms). The following high-order energy and dissipation functionals 
are fundamental to our analysis: 

S(t)=£(q,h)(t) 



■■=\ E iiA* 1/2 ^*iii5+5Ei(-^«) 1/2 «^- 1 » 8 - 26 ^ii 2 +5 E iim 1/2 (^+9^. 

a+26<5 " b=0 a+2&<6 

+ E ll(i-M) 1/2 ^||i s +i E II (i - /*) x/a + Saa** - «) ni s (i.ii) 

|<5|+26<5 |5|+2t<6 

and 

V(t)=V(q,h)(t) 

: = E llM^^iilj+EK-W^iJJ- 1 ^- 2 ^!!^ E ||M 1/2 (^ft+0°#M 

a+26<6 6=0 a+2fc<5 

+ E II(i-m) 1/2 ^M||j+ E ll(i-M) 1/2 (a 5 a t b gt +9 Q -a t b vi/ t . w )||^. (1.12) 

|a|+2h<6 |a|+2b<5 
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Note that the boundary norms of the height function are weighted by ^—d^q- We thus introduce 
the time-dependent function 

X (t):=m{(-d N q)(t,x)>0, 
which will be used to track the weighted behavior of h. We will show that 8 is indeed equivalent to 

3 3 

Ell^ll^(o)+xWEl 9 *^- 21 (r)' 

1=0 1=0 

and that T> is equivalent to 

2 2 

lkllH 8 .5 (n) +X]ll 9 ie*ll^-«(n)+xWEl 9 t +1/l l^- 2 '(r)- 

i=o i=o 
The elliptic operator in the parabolic equation (|1.9b ) for q has coefficients that depend on A = 
[D^]" 1 , which in turn depend on h; hence, the regularity of q is limited (and, in fact, determined) 
by the regularity of h on the boundary T. Since the regularity of h is given by norms which are 
weighted by the factor x(*)i a naive application of elliptic estimates would thus lead to the crude 
bound 

\\d l t q\\l 5 -2i<^), (1-13) 

which could a priori grow in time. However, by using the fact that lower-order norms of q have 
exponential decay (in time), estimate (|1.13[) can be improved to yield 

h\\l 5 + \\d l t q t \\t- 2l <e-"< t £ + V, 1 = 0,. ..,2 (1.14) 

for some positive constant 7>0. This is one of the essential ingredients of our analysis, as (11.141) 
will be used to control error terms arising from higher-order energy estimates in Section [3J 

In order to capture the exponential decay of the temperature q, we introduce the lower-order 
decay norms: 

2 1 2 

^(*):=^(EH^(*)ll^--W + Ell 9 Ml^--w)> 0(*) : =Ell#«(*)llH-»(n). (i-is) 

b=0 b=0 b=0 

with the constant /? denoting a positive real number given by 

f3~2X 1 -r h (1.16) 

where Ai is the smallest eigenvalue of the Dirichlet-Laplacian on il = i?i (0) and rj is a small positive 
constant related to the size of the initial data, which will be made precise below. Note that the 
smallness of Ep in particular implies an exponential decay (in time) estimate for the 7? 4 -norm of 
the temperature q(t). 

1.4.5. Taylor sign condition or non- degeneracy condition on qo- With respect to go = Po 0< ^i condi- 
tion (|1.2[) becomes ini x( z r [—dNq (x)} > A > on T. For initial temperature distributions that are not 
necessarily strictly positive in f2, this condition was shown to be necessary for local wcll-poscdness 
for (| 1.1 1) (see 29, SHIIH]). On the other hand, if we require strict positivity of our initial temperature 
functiorQ 

<7o>0 inO, (1.17) 

then the parabolic Hopf lemma (see, for example, [21]) guarantees that —djyq(t,x) > for < t < T 
on some a priori (possibly small) time interval, which, in turn, shows that £ and T> are norms for 
t>0, but uniformity may be lost as t— >0. To ensure a uniform lower-bound for —dNq{t) as t— >0, 
we impose the Taylor sign condition with the following lower-bounc0: 

-dNqo>C / q (fidx, (1-18) 
Jn 



2 Condition fll . 1 Tt is natural, since it determines the phase: Q(t) = {q(t) > 0} 
3 When h 

= (— sin9,cos0) 



^When ho ^ 0, the unit normal to the initial surface To is given by N = ^ +h " >s J> = J>Z where £ = (cos0,sin0) and 

v / (l+ho) 2 +9e'>o 



() 



MAHIR HADZIC AND STEVE SHKOLLER 



Here, ip\ is the positive first eigenfunction of the Dirichlct-Laplacian, and C > denotes a universal 
constant. The uniform lower-bound in ()1.18p thus ensures that our solutions are continuous in 
time; moreover, (|1.18[) allows us to establish a time-dependent optimal lower-bound for the quantity 
x(t) = 'm{ xe r(— djyq)(t,x) > for all time t>0, which will be crucial for our analysis. 

1.4.6. Compatibility conditions. The definition of our higher-order energy function £, restricted to 
time i = 0, requires an explanation of the time-derivates of q and h evaluated at t = 0. Specifically, 
the values qt\t=o, Qtt\t=0i n t\t=o an d h t t\t=o ar< 3 defined via space-derivatives using equations (|1.9a[) 
and ([OH]) . To ensure that the solution is continuously differcntiablc in time at t = we must impose 
compatibility conditions on the initial data (such conditions arc, of course, only necessary for regular 
initial data). By restricting the equation (|1.9ap to the boundary at time t = and using the fact 
that qt(0) — on T and that Af\t=o = 5f, where denotes the Kronecker delta which equals 1 if 
k = i and otherwise, we obtain the first-order compatibility condition 

Aq = (d N q ) 2 onT. (1.19) 

Upon differentiating (|1.9a[) with respect to time, and then restricting to T at t = and using (|1.19l) , 
we arrive at the second-order compatibility condition 

A 2 q = A\d N qo\ 2 +2d N (Aq -\d N q \ 2 )d N q -2\d NN q \ 2 onT, (1.20) 

where we have used that h t (t,0) =v [N(e)-T(9)h e (l + h)- 1 }. 

We note that our functional framework only requires specification of two higher-order compati- 
bility conditions (the condition <?o = on T being the zeroth-order condition). 

1.4.7. Main result. Our main result is a global-in-time stability theorem for solutions of the classical 
Stefan problem for surfaces which are nearly spherical and for temperature fields close to zero. The 
notion of "near" is measured by our energy norms as well as the dimensionless quantity 

K:=^. (1.21) 

IMIo 

as expressed in the following 

Theorem 1.2. Let (go,/io) satisfy the Taylor sign condition il.18]) . the strict positivity assump- 
tion ([7777]), and the compatibility conditions hi. 19)) . \1.20\) . Let K be defined as in \1.21\) . Then 
there exists an eo > and a monotonically increasing function F: (l,oo)— such that if 

£(qo,h )< 1 ^ ) , (1-22) 

then there exist unique solutions (q,h) to problem (1.9\) satisfying 

sup £{q{t),h{t))<Cel, 

0<t<oo 

for some universal constant C>0. Moreover, the temperature q(t)—>0 as t—+oo with bound 

ll9ll^ ( n)<Ce-^, 

where /3 = 2Ai~ O(eo) and Ai is the smallest eigenvalue of the Dirichlet-Laplacian on the unit disk. 
The moving boundary T(t) settles asymptotically to some nearby steady surface T and we have the 
uniform-in-time estimate 

sup \h-h \ 4 . 5 <y/e^ 

0<t<oo 

Remark 1.3. The increasing function F(K) given in hi .22)) has an explicit form. For generic 
constants C,C > 1 chosen in Sections [7J| and [4] below, 

F(K) :=max{8K 2ceK2 ,C w (lnK) w K 2oe ^}. (1.23) 
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Remark 1.4. The use of the constant K in our smallness assumption \1 .22)) allows us to determine a 
time T = Tk when the dynamics of the Stefan problem become strongly dominated by the projection of 
q onto the first eigenfunction ifi of the Dirichlet-Laplacian. Explicit knowledge of the K -dependence 
in the smallness assumption U.22\) permits the use of energy estimates to show that solutions exist 
in our energy space on the time-interval [Q,Tk]- For t>Tj(, certain error terms (that cannot be 
controlled by our energy and dissipation functions for large t ) become sign-definite with a good sign. 

1.5. A brief history of prior results on the Stefan problem. There is a large amount of 
literature on the classical one-phase Stefan problem. For an overview we refer the reader to Fried- 
man [23], Meirmanov (33] and Visintin [37]. First, weak solutions were defined by Kamenomost- 
skaya [3T], Friedman [22], and Ladyzhenskaya, Solonnikov k, Ural'ceva [37]. For the one- 
phase problem studied herein, a variational formulation was introduced by Friedman & Kinder- 
lehrer [2_3] , wherein additional regularity results for the free surface were obtained. Cafarelli [5] 
showed that in some space-time neighborhood of points xq on the free-boundary that have Lebesgue 
density, the boundary is C 1 in both space and time, and second derivatives of temperature are 
continuous up to the boundary. Under some regularity assumptions on the temperature, Lipschitz 
regularity of the free boundary was shown by Cafarelli [6]. In related work, Kinderlehrer & 
Nirenberg [321 [35] showed that the free boundary is analytic in space and of second Gevrey class 
in time, under the a priori assumption that the free boundary is C 1 with certain assumptions on the 
temperature function. In [Jj, Caffarelli & Friedman showed the continuity of the temperature 
in d dimensions. As for the two-phase classical Stefan problem, the continuity of the temperature 
in d dimensions for weak solutions was shown by Caffarelli & Evans [5]. 

Since the Stefan problem satisfies a maximum principle, its analysis is ideally suited to another 
type of weak solution called the viscosity solution. Regularity of viscosity solutions for the two- 
phase Stefan problem was established by Athanasopoulos, Caffarelli & Salsa in a series of 
seminal papers [3j [4] . Existence of viscosity solutions for the one-phase problem was established by 
Kim [32], and for the two- phase problem by Kim & Pozar [33]. A local- in-time regularity result 
was established by Choi & Kim [TT], where it was shown that initially Lipschitz free-boundaries 
become C 1 over a possibly smaller spatial region. For an exhaustive overview and introduction 
to the regularity theory of viscosity solutions we refer the reader to Caffarelli & Salsa [9]. 
In [36], KOCH showed by the use of von Miscs variables and harmonic analysis, that an priori C 1 
free-boundary in the two-phase problem becomes smooth. 

Local existence of classical solutions for the classical Stefan problem was established by Meir- 
manov (see [3j5] and references therein) and Hanzawa [30] ■ Meirmanov regularized the problem 
by adding artificial viscosity to (|l.lb[) and fixed the moving domain by switching to the so-called 
von Mises variables, obtaining solutions with less Sobolev-regularity than the initial data. Similarly, 
Hanzawa used Nash-Moser iteration to construct a local-in-time solution, but again, with derivative 
loss. A local-in-time existence result for the one-phase multi-dimensional Stefan problem was proved 
by Frolova & Solonnikov [26], using L p -type Sobolev spaces. For the two-phase Stefan problem, 
a local-in-time existence result for classical solutions was established by PRUSS, Saal, & Simonett 
|41) in the framework of L p -maximal regularity theory. 

In a related work, local existence for the two-dimensional two-phase Muskat problem (with varying 
viscosity and density) was proved by Cordoba, Cordoba & Gancedo [13] and in three dimensions 
in |14j . Their methods rely on a boundary-integral formulation for the Muskat problem, together 
with the Taylor sign condition. In a subsequent work [12j . various global existence results were 
established. An overview can be found in [TO] . 

As to the Stefan problem with surface tension (also known as the Stefan problem with Gibbs- 
Thomson correction), global weak solutions (without uniqueness) were given by Almgren & Wang, 
Luckhaus, and Roger [2] [38] [45]. In Friedman & Reitich [25] the authors considered the Stefan 
problem with small surface tension, i.e. with a<^i 1, whereby (|l.lc[) is replaced by u = ctk, k denoting 
mean curvature of the boundary. Local existence of classical solutions was studied by Radkevich 
[43] ; Escher, Pruss, & Simonett [20] proved a local existence and uniqueness result for classical 
solutions under a smallness assumption on the initial height function close to the reference flat 
boundary. Global existence close to flat hyper-surfaces was proved by Hadzic & Guo in [28], and 
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close to stationary spheres for the two-phase problem in Hadzic [27] and PRUSS, Simonett, & 
Zacher [42] . 

A global stability result for the two-phase classical Stefan problem in a smooth functional frame- 
work was also established by Meirmanov [35J for a specific (and somewhat restrictive) perturbation 
of a flat interface, wherein the initial geometry is a strip with imposed Dirichlet temperature con- 
ditions on the fixed top and bottom boundaries, allowing for only one equilibrium solution. A 
global existence result for smooth solutions was given by Daskalopoulos <fc Lee [19] under the 
log-concavity assumption on the initial temperature function, which in light of the level-set reformu- 
lation of the Stefan problem, requires convexity of the initial domain (a property that is preserved 
by the dynamics). 

In [25], we established the local-in-time existence, uniqueness, and regularity for the classical 
Stefan problem in L 2 Sobolev spaces, without derivative loss, using the functional framework given 
by (|1.11[) and (| 1 . 12|) . This framework is natural, and relies on the geometric control of the free- 
boundary, analogous to that used in the analysis of the free-boundary incompressible Euler equations 
in Coutand & Shkoller [E][T6]; the second- fundamental form is controlled by a a natural coercive 
quadratic form, generated from the inner-product of the tangential derivative of the cofactor matrix 
a, and the tangential derivative of the velocity of the moving boundary, and yields control of the 
norm / r (— dNq{t))\d k h\ 2 dx' for any fc>3. The Hopf lemma ensures positivity of —d^qif) and the 
Taylor sign condition on qo ensures a uniform lower-bound as t — > 0; on the other hand, —dNq(t) 
as t— >oo, and so an optimal lower-bound for (—djyq(t)) for large t is essential to establish a global 
existence and stability theory. 

1.6. Methodology and outline of the paper. Our present work builds on our new energy method 
for the Stefan problem that we developed in |29j . We obtain global and uniform control of the 
geometry of the free-boundary by controlling the weighted boundary-norm sup tg [ T ] || vxW^lk for 
all t > 0. We are thus able to track the location of the moving free-boundary and measure its 
deviation from the initial state; this geometric control is strongly coupled to, and dependent upon, 
the exponcntial-in-timc decay of the temperature function to zero. 

There exist infinitely many steady states for the classical Stefan problem: for any sufficiently 
smooth hypersurface rcIR d , the pair (p,r) = (0,r) forms an equilibrium solution of the Stefan 
problem (jl.lj) . This abundance of possible attractors for the long-time behavior of the solution T(t) 
creates a conceptual difficulty in approaching the question of "asymptotic" convergence. 

We address the temporal asymptotics by requiring our initial surface to be a small perturbation 
of the reference sphere. We use the energy spaces introduced in [25] ; moreover, we do not expect to 
observe any decay for the height of the moving surface in this norm. Rather, given the expectation 
that the solution does converge to some nearby shape (so that h remains small), we expect the 
temperature q(t) to converge to zero exponentially fast, since it is a solution of the nonlinear heat 
equation (|1.9a[) . Returning to the definition of the energy space £ given in (|1.11|) . wc immediately 
encounter a potential problem for global-in-time estimates; specifically, the coefficient —dNq(t) in 
the energy expression / r (— dNq(t)) \B 6 h\ 2 d9 is also expected to decay as t — >oo and it is a priori 
unclear how to uniformly-in-time control the regularity of the boundary height function h. To 
understand the relationship between the decay of q(t) and the smallness of £ , we will analyze the 
dynamics in three different and coupled regimes. 

High-order energy estimates. We do not expect the height function h(t) to decay to as t— >oo; 
rather, we expect h(t) to remain close to the initial height function h$. Assuming, without loss of 
generality, that ho = 0, to guarantee the smallness of h — ho = h we will prove that 



where P is some polynomial function of the low- norm Ep. The above estimate yields an a priori 
bound on £ if e, 8 and £ (0) are sufficiently small. 





(1.24) 
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However, to close the higher-order energy estimates and thus obtain (jl.24p . we must contend with 
a very problematic integral (or error term) given by 

Af :=- [ [ d N q t \d 6 h\ 2 d6dt. 
Jo Jr 

Driven by intuition from the linear heat equation, we expect dnqt to decay exactly as fast as —d^a. 
Comparing Af to the energy contribution J r (— 9jvq)|9 6 /i| 2 above, we note that Af cannot be controlled 
by £ , as it is the same order as £. Hence, to bound Af, we prove that after a sufficiently long time 
has elapsed, the quantity dNqt turns strictly positive and hence Af can be bounded from above by 
zero. In Lemma 14.21 we will quantify the meaning of "sufficiently long" time t = Tx from the previous 
sentence, expressing it as a function of the ratio K= 1 1 <?o 1 1 4 / 1 1 go 1 1 o - 

More precisely, we break the total time interval into a (possibly long) transient interval [0,T^] 
and [Tr-,oo). On the transient time-interval [0,7V] we do treat Af as an error term, and by choosing 
£ (0) sufficiently small, a straightforward application of a Gronwall-type inequality verifies that the 
interval of existence is greater than Tk, as explained in our proof of the main theorem (given 
Section 14.41) . The bound for Af grows exponentially with time, and as such, cannot be used to 
establish global-in-time estimates. Instead, a significantly more refined analysis is employed on the 
time- interval [Tr- ,oo), wherein we prove in Lemma 14.21 the negativity of Af for t = Tx and then use 
a maximum principle- type argument to guarantee the negativity for all t>Tx- 

Exponential decay-in-time of the temperature function q. The last inequality in (| 1 . 24|) holds only if 
Ep itself remains small; in fact we will prove that as t— s-oo, ||<z(t)||! has the nearly optimal decay 
rate 

e -(2A 1 -Ce )t ) ( L25 j 

where Ai denotes the smallest eigenvalue of the Dirichlet-Laplacian on the unit disk. Moreover, the 
parabolic estimate we prove, will be roughly of the form 

e -pt/2 

d t Ep + D<C(e + \\qo\\ 4 —r U ^)D 7 (1.26) 

where the norms Ep and D have been defined in (|1.15l) . A nice consequence of our analysis is that 
the potentially growing term, ^7^372 > hi f ac t remains small and decays in time. Next, we explain 
why this is true. 

Lower bound for the velocity of the free boundary. We may think of the presence of the denominator 
x ( t )i/2 ha the estimate (|1.26p as a possible obstruction to controlling the regularity of h and thus 
potentially preventing uniform ellipticity bounds for the parabolic operator (|1.9a[) . To deal with this 
issue, we need a quantitative lower bound on the decay rate of xW- Moreover, this lower bound 
has to favorably compare to the size of e _/34 . With some extra work, such a Hopf-type inequality is 
implied by a result of Oddson [3D] , which leads to the lower-bound 

x(t)>c ie - {Xl+ceo ^, (1.27) 

where c>0 denotes a generic constant, and as before c\ = J^qoPi is the first coefficient in the 
eigenfunction expansion of the initial datum qo with respect to the L 2 -orthonormal eigenbasis of 
the Dirichlet-Laplacian on the unit disk. Finally, combining (|1.25[) and (|1.27p . we will show in 
Lemma 12.31 that for small initial data, 

p -pt/2 

11*114^)172 ^v 7 ^ 4 . (1.28) 

for some positive constant 7* . 

Future work. In forthcoming work, we plan to address the Stefan problem on arbitrary domains 
diffcomorphic to the unit ball, as well as the case of the two-phase Stefan problem. Furthermore, we 
believe that the techniques developed in this paper should provide a general and robust framework 
for addressing the global stability questions for related free boundary problems in absence of surface 
tension. 
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Plan of the paper. In Section [2] wc introduce the bootstrap assumptions and obtain various a priori 
estimates, that allow us to control low norms of the boundary function h as well as the decaying 
low- norm Eg, and also establish the equivalence between the energies and the norms as mentioned 
earlier in the introduction (Section I2.5|) . In Section [3] we state energy identities and then perform 
the energy estimates. Finally in Section U we prove the main theorem. In Section 5, we discuss the 
modifications required for the analysis in three space dimensions, and for initial height functions 
ho 0. Appendix [A] is devoted to the proof of the energy identities stated in Section [3l The very 
short Appendix B provides a simple proof for the upper bound of dNqt- 

2. Bootstrap assumptions and a priori bounds 

Let us assume that the solution (q,h) to the Stefan problem (|1.9p exists on some time interval 
[0, T] , T > 0, which is guaranteed by [55] . With the positive constant eo < e <C 1 to be specified later, 
we make the following bootstrap assumptions: 

(smallness) sup £(s)+ [ V(s)ds<e 2 , sup Ep(s) + [ D{s)ds <CE a {0) , (2.29a) 

0<s<T Jo 0<s<T Jo 

(lower-bound) x0) £ c 1 e~ {Xl+1 ' /2)t , (2.29b) 

where we the definitions of £ , T>, Ep, and D are provided in p. lip . (|1.12[) . and (|1.15[) . respectively. 
With fi given in (|1.16p . /3 = 2Ai — n, the bootstrap assumption (|2.29b|) can be written as x(*) ~ 
cie~^' 2+ri ' t . Moreover, n > is a fixed small constant and it will be shown in the proof of the main 
theorem, Section \AA\ that 77 must be chosen smaller than 1 / v ChiK for some universal constant 
C. Note that since Ep(0)<e 2 , (|2.29ap implies the decay estimate \\q\\ 2 < e 2 e~@ t . Recall that the 
constant c\ in the estimate (|2.29bp is defined as f Q qo(x)ipi(x)dx. 

We now briefly explain the logic of the proof of global existence that will be carried out in 
Section @] If T is defined to be the maximal time at which the solution (q,h) exists and satisfies 
the bootstrap assumptions, the first objective is to show that the bootstrap assumptions (|2.29al) 
and (|2.29bp yield an improved smallness and lower-bound estimates at time T . If T were finite, by 
the local-in-time wcll-posedness theory and continuity of our norms we can extend the solution to 
an interval T+T*, while preserving the bootstrap assumptions (|2.29al) and (|2.29h>[) . thus arriving 
at contradiction to the definition of T. Hence T must be infinite. 

It remains to show that for e chosen small enough, the smallness and the lower-bound estimates can 
indeed be improved. In Corollary 12 . 141 we will show that the assumption (|2.29bp is in fact improved, 
and in Lemma [4.11 we show that the assumption on Ep + J Q D in (|2.29a|) is also improved. Finally, 

in Section [4.41 we will prove that the smallness of £ + J^V a ssumed in (|2.29a|) is also preserved. 
Thus the smallness regime introduced through (|2.29a|) — (|2.29b|) will be shown to remain preserved 
by the dynamics of (|1.9p for e > chosen sufficiently small. 

2.1. Poincare-type inequality. Because the first eigenfunction (p% of the Dirichlet-Laplacian is 
positive in fi, while the remaining eigenfunctions oscillate about zero, it will be necessary to introduce 
a constant into our estimates which gives a measure of the initial temperature distribution in the 
first mode of the dynamics. To this end, we will make use of the following 

Lemma 2.1. For fc>3, let f G H k (tt)P\HQ(tt), ,/:fi— s-R + be a strictly positive function on the 
interior of Q. Let ipi be the first eigenvector of the Dirichlet-Laplacian on the unit ball i?i(0) = O. 
Then there exists a universal constant C* such that 

\\f\\l<C*{ [ /(:r)^(:r)rf 3 ;)||/||3. 

Proof. We have that 

r ff x \ r 

/ f 2 dx< max — --^ / fifidx. 
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Since — ^^-(x)>c>0 for all x£T, the higher-order Hardy inequality (Lemma 1 in [T5]) together 
with the Sobolev embedding theorem shows that 

max l^)<C 



x£Q fi(x) 



<c\\f\\ 



which proves the lemma. 



□ 



Corollary 2.2. Let go € -ff 4 (il) CiHq (fl) with go > in f2. We consider the eigenfunction expansion 
Qa = '}2 r jLi c j l Pj °f 1o with respect to the L 2 -orthonormal basis {ipi,tp2,...} consisting of the Dirichlet- 
Laplacian eigenf unctions on the unit disk -Bi(O) = Q. Then, if jj*J^ < K , it follows in particular that 

M <K , j = l,2,... 

Cl 

Lemma 2.3. If the bootstrap assumptions \2.29a\) . \2. 29b)) hold, then 



x(t )l/2 

where 7 = 5 — r\ > 0. 
Proof. By (|2.29b|) . we have that 
Epit^e-^ 2 



x (i)i/2 



(2.30) 



-Pt/2 



x{t) l/2 



< C 



1/2 



e -(Ai/2+7)/4)t 1/2 



< c e -^ 2 « < CK\\q \\l /2 e-^ 2 < CV~ee~ 
1/2 



7*/2 



where we have used the fact that c^ 2 > -^72 Ikollo^ anc ^ IkolU^S-^lkollo- We have also used the 
bound if || go II 4^ 2 < Cy/e (since e < e), as well as the smallness assumption (|1.22[) so that if || go II 4^ ^ 
Keo/FiK) 1 / 2 <Ce. Note that 7 is explicitly given by 7 = (§ -77) > 0, and that t?<Ai/2. □ 

2.2. A priori bounds on h. 

Lemma 2.4 (Suboptimal decay bound for ht). Under the bootstrap assumptions \2.29a\) and \2.29b\) . 
the following decay bound holds: 

\h t \2. 5 <ee^ 1 / 2 . (2.31) 

Proof. Differentiating equation (|1.5[) . the Sobolev embedding theorem together with the fact that 
h>0 (by the maximum principle) show that 

Ni < Mwi.~ + N2M2 + N1M1N1 

<\v\2 + \h\2\v\ 2 + \h\i\v\i(\h \i+t sup \h t \i), 

0<s<t 

where we have used the fundamental theorem of calculus for the last inequality. Using the bootstrap 
assumption (|2.29a[) . we see that \v(t)\ < e - ' 3 *, while thanks to Lemma [2.31 and the fact that y/£< 
£o<e, 



\h\2\v\ 2 <Vx\h\ 2 



< 



'E. 



Hi 



lV£e-^l 2 <ee-^' 2 



Hence, 



sup \h t \ 1 <ee-^ 2 + ee^ t / 2 (l+ sup |^|i), 

0<s<t 0<s<t 



and with e > sufficiently small, we see that 



sup Ni<ee- 7 * /2 <e. 

0<s<* 

Taking more derivatives of (11.51) . the Sobolev embedding theorem shows that for A: = 2, 3, 



he 



he 



(2.32) 
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where we have again used the fact that h > 0. Since 

I he 

for some polynomial function P. and since \h\k < |/io|fc + £sup 0<s<t \h t \k, we see that 



<\h\ k+1 {l + P{\h\k-x)) t fc = 2,3, 



|-f-| k <\h\ k+1 (l + P(t)P( sup 1/itU-x)). (2.34) 

We now use (j2~M]) and ([2331) to infer that 

\ht\k<\v\ k {l+ sup |/H| 2 ) + |/i|fc+iNi(l + P(t)P( sup Ifctla)), (2.35) 

0<s<t 0<s<t 

where we have used \hg\i ;$£sup 0<s<t \h t \ 2 . Interpolating between k = 2 and k = i yields 

|ftt|2.5<|«|2.5(l+ sup \h t \ 2 )+\h\2. 5 \v\i(l + P{t)P{ sup \h t \ 2 )). (2.36) 

0<s<* 0<s<t 

and as above, Lemma [2.31 provides us with the inequality | /i| 2.5 |f 1 1 ee -7 '/ 2 , which together with 
the bootstrap assumption (|2.29a[) shows that 

sup \h t \ 2 . 5 <ee-W 2 (l+ sup \h t \ 2 . 5 )+ee~^ 2 {l + P(t.)P( sup \ht\2)) 

0<s<t 0<s<t 0<s<t 

and therefore with e > sufficiently small, 

sup N 2 . 5 <ee-^/ 2 + e e - 7t / 2 (l + P( sup \h t \ 2 )), (2.37) 

0<a<t 0<s<t 

where the polynomial P[t) has been absorbed in some universal constant due to the exponentially 
decaying factor e -7 '/ 2 . On the other hand, the inequality (|2.35[) with k = 2 together with the estimate 
(|2.32[) shows that \h t \ 2 <e so that with (|2.37[) , we conclude the proof. □ 

Remark 2.5. Note that the estimate \2.31\) can be stated more precisely, by keeping track of constant 
c\ on the right-hand side, in which case, 

Ihth.^^V^e-^ 2 . (2.38) 

The proof follows from the last line of the proof of Lemma \2.4\ since ^(O) 1 / 2 < K 2 c\, due to the 
bound \\q\\ 4 <K\\q \\ <K 2 c 1 . Note that y/1 on the right-hand side of \2. 30\) can be replaced by ^fc{ 
for the same reason. 

Lemma 2.6 (Smallness of the height function). Let c\~ J^qoipidx and suppose that the bootstrap 
assumptions V2.29a}) . 12. 29b)) hold. For e > taken sufficiently small, 

sup \h(s)\ 4 . 5 <V~t, (2.39) 

0<s<t 



while for lower-order norms, 



sup \h(s)\ 2 . 5 < Cl and sup \h(s)\ 4 < e 1/2 cj /4 . (2.40) 

0<s<t 0<s<t 

Proof. Observe that 

\h\ 2 2 5 <2 f \h\ 2 , 5 \h s \ 2 , 5 ds< sup \h(s)\ 2 , 5 f \h s \ 2 . 5 ds< sup \h(s)\ 2 . 5 f e^^Te" 7 '/ 2 , 

JO 0<s<i JO 0<s<t Jo 

where we have used (|2.38[) in the last bound. Taking the supremum over the time interval [0,t] we 
deduce 

sup |^(s)| 2 .5<e 1/2 y5T. 

0<s<t 

Using the well-known interpolation estimate (see, for example, pQ) 

|/k.<|/|f|/|^ 9 , e = ^,l<k<m, (2.41) 
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with fc = 3, 1 = 2.5, m = 4:, and the fact that \s/xd ht\a is bounded by £ , we have that 

where 7* = ~ 57+£(§ + §) = ~ + if > ^ (by definition, 7 = — | — 77). As a consequence, 
Wl< / |ft|sN3< sup |ft(s)| 3 [ \h t (s)\ 3 ds<e sup |/i(s)| 3 . 

JO 0<s<t Jo 0<s<t 

Upon taking the supremum over the inetrval [0,t], we finally have that 

sup \h(s)\ 3 <e. (2.42) 

0<s<t 

We can now improve the decay result of Lemma \2~M first for the quantity \ht\2- Simply using the 
bound (|2.42[) . exactly as in the proof of Lemma 12.41 . we infer the improved estimate 

N 2 <|M| 2 . 5 (1 + N3)<cie-^ 2 . (2.43) 

As an immediate consequence, we obtain the smallness bound for sup 0<s<t | ft.(s) [4: 



I 



d 4 h\ 2 d6 = d 4 hd 4 h t d6ds = d 6 hd 2 h t d9ds 
jo Jr Jo Jr 

< /'|a 6 /i|o|«9 2 ^|ods< [\4^Lcie- Ps/2 )d S 
Jo Jo X\ s ) 1 

< f eV3Te" 7s/2 rfs<eVcT- 
Jo 



(2.44) 



Note that (|2.44j) . in particular, implies the second bound in (|2.40[) . Next, we establish the a priori 
smallness of sup 0<s<t |/i(s)|4.5. Thanks to (J2T44J), we improve the decay bound for | /z 4 1 2.5 in an 
analogous fashion to the improved decay estimate (|2.43j) for \h t \2- We obtain | /i t 1 2 .5 ^cie _,3t / 2 . The 
first bound in (12.401) now follows from the fundamental theorem of calculus and the previous bound. 
A straightforward interpolation argument for fractional Sobolev spaces on the unit circle V, shows 



\h\i 6 < / \h\ 6 \h t \ 3 ds. (2.45) 
Jo 

Using the interpolation estimate f|2 .41 [) . with 1 = 2.5, k = 3, and m = 5, we see that 

\h t \ 3 <C\h t \^\h t \l /5 . (2.46) 
Using (|2.46[) with (|2.45[) . and using the above bound on | /i. t 1 2. 5 ? yields 

\h\h< [\h\ 6 \h t \^\h t \l/ 5 ds< f*!L#\-W*\h$*<U 
Jo jo X\ s ) ' 

<ec\' W fe-^h^ds, 
Jo 

where we also used the bootstrap assumption (|2.29b[) . One checks that 7 = — ^ + {4 + 2) = lb~P ~ 2 > 
0. We thus have 

I ds. 



l^,<^ /1 7 o e-^ 2 x(e-^ 2 |^r) l 
Holder's inequality with p = ^ and q = 10 then shows that 

m, < *cr °({( e --/ 2 ) io/9 ds ) 9/10 (/(e-^i, t i 2 ds ) l/l0 
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where the last inequality follows from the definition of 7 above, the bootstrap assumptions (|2.29bt) 
and (|2.29a|) . and the estimate 

t r t , 

e-^ s \h t \ 2 5 ds< / —e-^ s+ ^/ 4+ ^/ 2 >mi{-d N q(s))\h t \lds 



c i 
t 



L e -WW)'w£(-d N q(8))\ht\lda 
ci r 

<-/ f(-d N q{s))\d 5 h t \ 2 deds<-. 



< 



ci J Q J r ci 

□ 

2.3. Differentiation rules for A. Since A = [D^] -1 , it follows that 

d t A\ = -A k r w r >s At ; M\ = - A k r 8* r , s A\ \ 

In particular, a simple application of the above identities and the product rule imply that for any 
given a, be N, 

B a d b t A k i = -A k r B a d b t ^ r , s At + {B a d b , A k } , (2.47a) 
{<9 m <9",Af}:= ]T au,d l d l i(A k Al)d m - l d?~ l '-V r , s , (2.47b) 

l+V>l 

where the term {■, •} is the commutator error . Here the constants ai^y are some universal constants, 
depending only on to, n, I and V (where < / < to, < V < n). 



2.4. Estimates for V^ — Id and A — Id. Under assumption (|2.29aj) . the elliptic estimate (jl.7 
shows that on the time-interval [0,T], 



||V*-Id|U-(B l) <C , ||V*-Id||i.6<C|A| a (2.48) 

and for < s < 3, 

\\DH\\ s <C\h\ s +i. 5 . 

Estimate (|2.48|l implies that 

||^-Id|| L oo (Bl) = ||(Id-Vf)A|| L oo (Bl) <C||A|| L oo (Bl) |/ l | 2 ; 
thus under assumption (|2.29al) . 

\\A-U\\ L ~ (Bl) <C\h\ 2 (2.49) 
Note that (|2~48]l and (pT49j) together imply that for < s < 3, 

\\DA\\ s <C\h\ s+1 . 5 . 
Thus, with Lemma 12.61 we have proven the following 

Lemma 2.7. With the bootstrap assumptions i2.29a)) . 12. 29b\) andfore>0 taken sufficiently small, 

||Vtf-Id|| 4 + P-Id||4<v^- 

2.5. High-order derivatives of q. Because our energy function £(t) is formed using only tangential 
derivatives in space, the purpose of this section is show that radial derivatives of the temperature q 
are also bounded, and thus the full Sobolev norms of the temperature q are controlled by our energy 
function, as was explained in the introduction. 

We will make use of the heat equation and its time-differentiated variants: 

ft-A*? = / 0) (2.50a) 
qtt-^qt = h (2.50b) 
qttt-A-9qtt = f2, (2.50c) 
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where A*=A|g|- \A k and where the forcing functions /o,/i,/2 are given by 
fo = -V t -v, 

h = -{^-v) t +4{d t A k i q, k )„■ +9 t A? (4 3,* ),i > 

/a = - (*t ■ «)« + 2^ (3, A* ft )fc )„■ +29 t ^ (A k g t , fc +25 t A? ($i4* g, fc ) „■ 

+ ^'(^ fc g t ,j k ) JJ -+^(^?, fc ), 3 -. 
We will repeatedly make use of the following elliptic estimate: 

Lemma 2.8 (Elliptic regularity with Sobolev-class coefficients). Let q denote the unique Hq(Q.) 
solution to 

—A^q = F in tt, 
q = on dfl . 

Suppose thatk>l, FeH k - 1 (Q), and AeH k (n) satisfying A^O- > A|£| 2 f or a^£eR 2 for some 
A>0. Then 

\W\\H k + 1 (n)<C \\ F \\ H k - 1 (n) + \\M P Hk(n)\\ F \\L"(n) 

for some power p>l. 



Proof. We provide the details in the course of the proof of Lemma 12.91 □ 

Lemma 2.9 (Bounding d\q, Z = 0,l,2,3, by £(t)). With the bootstrap assumptions ¥2. 29a)) 
and 12. 29b)) holding, and with e > sufficiently small, there exists a constant C* such that 

lltello + lk«ll2 + ll*lll+lkll 2 <^. 

Proof. Stepl. Estimating \h U t\o.5- We denote by X(t) the quantity H^tttjlo + Win -111 + \\ c li 111 + Iklli- 
Twice time-differentiating (|1.9d|) . we find that 

h m = V tt -N-\^-}„V-T-2\^-\ t V t -T-^-V tt -T. (2.51) 

l l + h ht l l + h h l + h v ' 

By the normal trace theorem (see, for example, equation (6.1) in |18j ). 

|««-JV|o.6<||9wtt||o + l|div«tt||§. 

Note that 

divwt t = (div* v) tt + ((div- div*)w) tt = {q t + v- %) tt + [(A* - S k )v\ k ] tt = q ttt + * ttt -v + K, (2.52) 

where the remainder 1Z reads 

K = 2*« -v t + * t - v tt + (A k - S k ) tt v\ k +2(A k - 5 k ) t «* fc + (A k - 5 k )v^ k . 

From Lemma [2.61 and [ 2.31 we obtain the estimate \\7Z\\q < e£ + eX. Thus, returning to (|2.52[) and 
using that \\qttt + ^ttt • v|lo < £ by (|l.lll) . we get ||divu tt ||Q<£ + eAf and consequently 

\vtt-N\o. 5 <£ + eX. (2.53) 

As for the last term on the right-hand side of (|2.5ip . we use the tangential trace theorem (see, for 
example, equation (6.2) in [TB]) to infer that 

l«« • t| ^ ll^wtt Ho + ||curlw«||g. 

Since curl$u = (recall v = — Vpo^/), we have curlw tt = [(curl — curhj )v] t t- By a similar inequality 
as above, using Lemmas 12.61 and l2.3[ we obtain the bound || [(curl — cmly)v]u\\o < c£ + eX . Together 
with (f2~53| and \\dv u \\l < £ , this leads to 

\v t fT\ .5<£ + eX. 

Together with the smallness of hg and hg t from Lemma l2~6l the bound \^fxdh tt \ \ < £ and Lemma [2~3l 
we finally infer from (|2.5ip that 

\h m \o.5<£ + eX. (2.54) 
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Step 2: L 2 estimates for d\q. By the triangle inequality and the definition of £(t), we have 

that for ; = 1,2,3, 

\\dlq\\t<\\d l tq + d l t *-v\\t + \\d l t *-v\\ 2 
<£(t) + \\dl*-v\\l 

<S(t) + \\v\\l\\dlnl<m + e 2 \h m \l 5 
<£{t) + eX, 

where we used the Sobolev embedding theorem and (|2.54[) . 

Step 3: H 2 estimate for qtt- We consider the elliptic equation — A^g = /o — q t . We note that 
Lemma 12771 ensures that A 1 ] Wiii , > h\£\ 2 forall£eR 2 . Given that ||/ — q t \\Q <£ , elliptic estimates 
show that This, in turn, implies that — ?tt||o~£j an d elliptic estimates then show 

that Hftlli^^- Hence, we have that H/2 — QtttWo ^£ + eX, and once again use elliptic estimates to 
conclude that \\qttW2 <£ + e%- 

Step 4'- H A estimate for q t . Since ||/o — qtWl 1^£, Lemma [2~8l shows that Hoill^f; thus, — <Ztt|||^$ 
£ + eX. Another application of Lemma l2~8l together with Lemma l2~7l thcn shows that \\qt\\ 2 <£ + eX. 

Step 5: H 6 estimate for q. The elliptic estimates in Steps 3 and 4 made use of Lemma l2~7l To 
obtain the H 6 estimate for q requires us to improve the elliptic estimate in Lemma 12.81 to be linear 
in || Vx*l|6- To this end, we write A jk =AjA k and rewrite (|2.50a ) as 



- (A ]k q, k ) ,j = -qt + fo — A\ ,j A k q, 



k ■ 



(2.55) 



Letting d a act on (|2.55j) . we find that d a q satisfies 
-'■A ij (d a qh 



8 a (*fv + q t )+ J2 C a p[(d^A ij ){d a -^q), 

0<f3<a 

Y C aP ¥(AUA k )d a -?q, 3 , 



0<f3<a 

where C a p are constants from the product rule. Multiplying this equation with d a q and integrating- 
by-parts, using the fact that B a q = on dil and that A> 1/2, we find that 

0<f3<a 

+ E C a0 \\d^Al 3 A k )d a ^q, 3 \\ o \\d a q\\o + 

O<0<a 



l\\d a q\\ 2 1 <\\d a - 1 ^ t -v + q t )\\ \\d a+1 q\\o+ ]T C a0 {B A ij )(B a ~ q), 



\B a q 



i 



^, 3 A k )Bq 



l-J 







\B a q\ 



(2.56) 



Let us examine the second term on the right-hand side of (|2.56[) . By Young's inequality, for 5>0 

2 





Y C a p (B^A ij )(B a ~ q),j \\8 a q*\\o<S\\d a q\\ a 1 + C 6 Y C ^ d?AB a -PDq 

0</3<a 0</3<a 

where C S = C/S. Since B 5 A^B 5 D^P(A) + B 4 D^P(dD^,A) + B 3 D^P(B 2 D^,BD^ : A) : it thus fol- 
lows that for a = 4 or 5, 



\\d a A\\ <C\\d a - 2 (^ -e)\\i<C 



,\Vxh\h ,£ 



< . 



(2.57) 



X X 

The linear inequality (|2.57[) shows that our bootstrap assumptions (|2.29a[) and (|2.29b[) imply that 
the map In— > A is linear with respect to these high norms. 

We first consider the case that a = 4. Prom (|2.57|) when a = f3 = 4 



I d a A Dq\\ 2 <-^- £<ee-^£. 



X 



(2.58) 



The Cauchy-Schwarz inequality, together with the Sobolev embedding theorem, shows that 
W&A dDq\\ 2 has the same bound. Next, \\d 2 A B 2 Dq\\ 2 + \\dA d 3 Dq\\ 2 < \\V\\l\\q\\l < ee~^£ < ee~^£. 
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The first, third, and fourth terms on the right-hand side of (|2.56p are estimated in a similar 
fashion, so we do not provide the details. Hence, by choosing S > sufficiently small and employing 
Young's inequality, wc find that 



\<1 



l4 + EH^IIi<£ + ^. 



a<4 



To estimate radial derivatives, we use polar coordinates for the disc (with the usual basis e r and 
eg). Expressing the components of the matrix A as 



A- 

we may write 



A 6r A 6b 



&iv{A Vg) = r- 1 (rA rr q r ) r + r" 1 (.4 r %) r + ^{A^q^e + r' 1 {r' 1 A ee q e )e ■ 
It follows that 

-A rr d a q rr = r-\rAnrd a q r + r-\A re d a q r ) r + r- 1 {A re d a q r )g + r- 1 (r- 1 A ee d a q e )e 

~d a (^fV + q t )+ £ C^d^A^B^q)^ 

0</3<q 

- E CapBPfaiAfyff*-^, (2.59) 

0</3<a 

Let uj = {x£il : i<|a;|<l}. For a < 3, every term on the right-hand side has I/ 2 (w)-norm 
bounded by a constant multiple of £. Hence, it follows that 

Allowing -j^ to act on (|2.59|) . as many as three times, we conclude that 

IMIL<£ + ^. (2.60) 

We return to the inequality (|2.56p and consider the case that a = 5. Once again, we focus on 
the second term on the right-hand side, the first and third terms being similar (and easier). From 
(12.581) ||9 5 ^4 Dq\\g <ee~ 7 *£ . The Cauchy-Schwarz inequality, together with the Sobolev embedding 
theorem, shows that \\B A A dDq\\l + \\d 3 A B 2 Dq\\l + \\B 2 A B 3 Dq\\l < ee'^S. Finally, using (l2~60| . 
we conclude \\dA d 4 Dq\\l < e\\v\\l < ee _7 *£ . We conclude that 

\\ q \\l+J2\\9 a <i\\i<£+^. 

Then setting a = and letting act on (|2.59[) shows that indeed 

MIL, 



By using a smooth cut-off function whose support contains f2\w, we easily obtain the interior 
estimates, and find that \\q\\% < £ = eX. Recalling the definition of X and the estimates from Steps 
2, 3, and 4, we finally infer X<£, which concludes the proof of the lemma. □ 

Lemma 2.10 (Bounding d l t q, Z = 0,l,2,3 by T>{t)). With the bootstrap assumptions (2. 29a)) 
and H2.29b)) . and for e > sufficiently small, there exists a 7 > such that 

2 

Ell^H^ + lkllL^ee-^ + O. (2.61) 
1=0 

Corollary 2.11. With the bootstrap assumptions 112.29a)) . \2. 29b)) and a sufficiently small e>0, 

M\l. 5 + \ht\l<ee-^£+V 
with 7 = /3/2 — 7] as defined in Lemma[ 
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Proof of Corollary \2.11\ We write (jl.9b|) as 

v = Dq-(U-A)-Dq. 
Using the basic estimate from Lemma \2.7l we see that 

Il«ll5<(l + v^)lkl|6 + 4e-^ 2 ||*- e |j 6 , 
lkl| 6 <(l + ^)||9||7 + i?|e-^/ 2 ||*- e || 7 , 

so that an application of linear interpolation (see, for example, Theorem 7.17 in Adams fjQ) provides 
the inequality 

\\v\\l 5 <(l + V~e)\\q\\l 5 + E l3 \^-e\\l 5 . 
Using Lemmas 12.31 and 12.101 it follows that 

\Hh<(i+V$h\\h+-^ ) — x(*)ll*-e||l.B 



<(i+v^)lkllL+ PK £ t) m 

Next, using the formula (|1.5[) . we see that 

Ml<x{t)\h\l^ + e\v\\ 

which once again, thanks to Lemmas 12.31 and 12.101 is bounded by a constant multiple of e 2 e~ 7 *£ + 
V. □ 

Proof of Lemma \2.1U\ Step 1: H 1 estimates for d\q. We make use of the identity X7q = v-'V'$>. It 
follows that 

Vq tt = v tt ■ V* + 2v t ■ W t + v ■ V* tt 
Vqttt = vttt • W + 3v tt ■ V* t + 3v t ■ W tt + v ■ V* t « . 
Employing Holder's inequality and the Sobolcv embedding theorem, 

llVfell^KJ 2 , \h\t + \\vu\\ 2 N5 + ^ \Vxhu\l + ^ \Vxhut\ 2 Q <V 

where we have used Lemma 12.301 for the last inequality. We have similar estimates for q u , qt, and q 
so that 



£iiM<^ ( 2 - 62 ) 



1=0 

Step 2. H 3 estimate for q u ■ Just as in the proof of Corollary 12.111 we see that as a consequence of 
Lemma 12.91 

2 

J2\\d l t v\\ 5 - 2l <£. (2.63) 
;=o 

Returning to the equation (|2.50a ). we estimate —^t'V — qt in iJ 1 (f2). By the Sobolev embedding 
theorem together with Lemmas 12.41 and 12.31 

!l*t!l^-<ll* t ||3<v^e-^ 2 , (2.64) 

so that together with (|2~63| . ||* t • v\\\ < ee 7 *£ . Then, with with (f2J62]) . 

hWi^^-^S + V. (2.65) 
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Next, we return to (p.SOb 'l and estimate fi — qtt in By Lemma |2"H1 \\ty t ■ v t \\l < ee _7 *£ , 

while W^tt-vWl^^^Ep^ee-^S. The estimates flgjgg) and (|2~65|) then show that ||/i - q tt \\f < 
ee'^S+V so that 

\\qt\\l<ee-*£ + V. 

A similar estimate then shows that — <?ttt|li^ ee ~ 7 *£ + £> so that from ()2.50b ). 

hu\\l<ee-^e+V. 

Step 3. H 5 estimate for q t . From (j2~63f and (f2J4|) . we see that ||* t • u||§ < ee~ 7t £ +T>, so that with 
Lemmas 12.71 and 12.81 we have that 

ll<7ll5<ee- 7t £ + 2?. 

This, in turn, ensures that ||/i — g 4t 1 1 § < ee _7 '£ + V so that 

\\qt\\l<ee-^E + V. 

Step 4- H 6 - 5 estimate for q. We first look at the estimate (|2.56p with a = 5. We find that 
\\d 5 q\\l<\\*fv\\ 4 + \\q t \\ 4 + 

0</3<5 

+ E |^(4,i4 fc )^-^,J +|^(4,i^)^IL (2 - 66) 

0</3<5 

For the first term on the right-hand side, we note that with the Sobolev embedding theorem and 
Lemma 12.31 



ll*t^U<llv^lU^ + ll^l! 3 |MU 

<V~ee-> t/2 (\\^t\\k + \Mk) fc = 4,5. 
Using the estimate (|2.57p . we see that 

\\(B^A^)(d^q)J <V^ e -^ /2 (||Vx(*- e )l|6 + lkl|5) 

0<P<5 

The last two term on the right-hand side of (|2.66p are estimated in the same way so that 

l|5 5 9||i<^ e - 7t/2 (||Vx*t||4 + |IVx(*-e)|| 6 + ||«||4 + |k t || 4 ). 
Using the formula (|2.59[) . we find that 

lkl| 6 <^ 7 * /2 (||^j4 + |IVx(*-e)|| 6 + IK'||4 + !k t ||4) + lk t ||4. (2.67) 
The identical procedure with a = 6 then yields 

!kll7<^e- 7 * /2 (|IVx*tll5 + ll^(*-e)||7 + |K'||5 + lkj5) + lk t || 5 . (2.68) 
Linear interpolation between (|2.67[) and (|2 . 68[) , we have that 

ll^lle.s < V^e~^*/ 2 (|1 V^^*ll4.5-U || — e)|| 6 . 5 -K ||x,|| 4 . 5 + ||Q t || 4 . 5 ) + ||cz t || 4 . 5 



□ 
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2.6. Lower bound on x(t). The heat equation (|1.9a|) for q can be rewritten as 

qt-akjq,kj-hq,k = in (2.69a) 
g = 0on T, (2.69b) 
q(0,-) = <?o>0 in (2.69c) 
where the coefficient matrix a = (afcj)fcj=i.2, and the vector b= (61,62) are explicitly given by: 

a kj :=A*Ai; b k := A^M + A<?¥ t . (2.70) 

Wc first quote a theorem from |40j , that will play an important role in producing quantitative bounds 
from below for x(0- 

Lemma 2.12 (Oddson's Theorem 2 in [40]). Let q£ C 1,2 {Vl) be a supersolution to 12. 69)) in the unit 
disc f2 = i?i(0), and let <a< \ be the normalized ellipticity constant satisfying 

ajktijZk > a(au +a 22 ) |£| 2 
for any real vector ^=(^1,^2)- Moreover, let us introduce the quantities 

ko(T):= inf , f3{T) := sup b-x. 

ox[o.t] an+a 2 2 nx[o,T] 

Let denote the Bessel function of the first kind of order p and £0 its first positive zero. If we 
define 

la, kq 
then there exists a positive constant m satisfying 

q{t,x) > mpe~ xt , 

in -Bi(O) x [<t,oo[, where p stands for the distance from x to the boundary T and a is an arbitrary 
small time. 

Remark 2.13 (Optimal decay rate for solutions of the heat equation). If we set A = Id, then 
problem 12. 69)) turns into the initial-boundary value problem for the linear heat equation. In this 
case k = i } a = 1/2, f3 = 0, p = —1 = 0, and A = £o? where £0 stands for the first positive zero of 
Jo(£). In particular, if q heat denotes the associated solution, then the above lemma implies that 

Xhea t(t) := inf (-d N q heat (t,x)) > e~^\ 

which is the optimal decay rate in the case of the linear heat equation, as the lowest positive eigenvalue 
of the Dirichlet-Laplacian on the two-dimensional disk corresponds exactly to 

Corollary 2.14 (Lower-bound for x(t))- Under the bootstrap assumptions 12. 29a)) and 12. 29b)) with 
e small enough, there exists a universal constant C > such that 

x(0>cie- (Al+ ^ (t))t , 

where c\ = J^qoipidx is the first coefficient in the eigenfunction expansion of the initial datum qo with 
respect to the L 2 ortho-normal basis {(px,tp2,.-.} of the eigenvectors of the Dirichlet-Laplacian on 

Bi(0), i.e qo = ciipi+C2¥>2 + Moreover, A(i)>0 satisfying A(t)<Ce for some positive constant 

C. In particular, with e>0 sufficiently small so that Ce<i]/A, we obtain the improvement of the 
bootstrap bound 12. 29b)) given 

Proof. The proof of Oddson's Theorem 2 in [JUI (Lemma I2.12[) relics on the construction of a 
comparison function of the form v(t,r) = r~ fJ, J^(t;or)e~ :Kt , where A,/i,£o are given in the statement 
of Lemma [2.121 J M is a Bessel function of the first kind and r=|x| is the radial coordinate. The 
first property of v which is important for the proof is that v vanishes at the spatial boundary T and 
approaches it like c(l — r)e~ xt as r — > 1. This is a consequence of the fact that lim r ->i ^(l^r) = c ^ or 
some constant c>0, a well known property of Bessel functions. The second important property is 
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that v is a subsolution for (|2.69p (and it is constructed with the help of maximal Pucci operators as 
explained in detail in [40]). 

The goal is to prove that for any arbitrarily small time a > there exists a strictly positive 
constant 6(a) >0 such that q — 5v is a positive supcrsolution to the parabolic problem (|2.69l) on the 
time interval [cr,oo[. The desired lower bound for q then follows from the weak maximum principle. 

Since v is a subsolution, it follows that for any <5>0, q — Sv is a supcrsolution. The positivity 
of q — Sv at t = a follows from the parabolic Hopf lemma, from which we infer the existence of a 
constant 8(a) such that ^>S(a) uniformly over O. Note that we have used the fact that v(a,r) 
behaves like C(l — r) near the boundary T for some positive constant C. Therefore it follows that 
the constant m in the statement of Lemma 12.121 a priori depends on the time a > 0, and moreover, 
m is proportional to the lower bound for —dq/dN\ t=a on T. 

From the proof of the parabolic Hopf lemma (see for instance Theorem 3.14 in [21]), the value 
—dq/dN\t =a is proportional to the minimal value of the temperature q on a space-time region of the 
form K a := B\~ca x [ f7 /2,3cr/2], divided by a (which is roughly the distance of K a from the parabolic 
boundary of f2 x [0,2er]). Note that, unlike the elliptic case, we are forced to take into account the 
time-dependence of the solution and in particular the region K a cannot be chosen uniformly for all 
times, but only for times greater or equal some arbitrarily small a > 0. However, our solution is 
continuous all the way to t = and we do nevertheless obtain a lower bound for all times due to the 
Taylor sign condition; namely, due to (|1.18[) , 

—dNqo ^ 
-d N q = ci >ci. 

Cl 

Note however that if we define the dimensionless quantity L= {—dj^q^/ci >0 and assume no uni- 
versal bound on L from below, the only modification in the statement of the main theorem will be 
that the smallness condition on initial data (| 1 . 22|) will additionally depend on L. 

As to the bound on A, note that the exponent X — X((aij),(bi)) depends on the coefficients 
( a u)ij=i,2 an d (&i)j=i,2 through the relationship X — a^/kQ. Since kg and £o vary continuously 
as the coefficients are varied, it proves that A depends continuously on the coefficients a,ij,bi of 
the parabolic operator. On the other hand, by Remark 12.131 it follows A| Qi =5^-^=0 = Ai. As a 
consequence 

|A(0| = |A(t)-Ai|<C(||A-Id||^,||6||^) = O(||^ 2 (*-c)||^,||vI/ t || LM ). 

□ 



3. Energy identity and the higher-order energy estimate 

3.1. The energy identity. Much of our analysis is founded on basic higher-order energy identities 
for the classical Stefan problem. These identities provide the geometrical control of the evolving 
phase boundary, which in turn controls the decay of the temperature function; moreover, these 
identities explain our definition of the higher-order energy function £ and the dissipation function 
V. 

Proposition 3.1 (Energy identity). With R=l + h and Rj = i?J _1 , sufficiently smooth solutions 
to the classical Stefan problem satisfy 

where the error terms IZj, TZj, Sj, Sj, Qj, and Hj are given by \A. 1$ , KA.15\) , 
fA~M) . (A~¥fy , /COffp . and (A~2$ , respectively. 

The proof is provided in Appendix [3] 
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Remark 3.2. On the right-hand side of \3. 71% we have isolated the error term 

foopf = ^(-d N q t )R 2 j\B 6 - 2j d 3 t h\ 2 dx', 



(3.72) 



from the other boundary-integral error terms Qj and T-Lj; indeed, Gnopf can only be thought of as an 
"error term" on a transient time-interval, for after a sufficiently large time, we will no longer be 
able to control Gnopf u *a energy methods, and instead, we have to rely upon a Hopf-type argument 
to prove that Gu Q pf < 0. 

3.2. Energy estimates. To control some of the highest-order error terms in our energy estimates, 
we shall make use of the following technical lemma, whose proof is given in [15] and |16| . 

Lemma 3.3. Let Hz(£l)' denote the dual space of H^(Q). There exists a positive constant C such 
that 

WdF\\ i{ny <C\\F\\ Hi{n) forFeHHn). 

As a consequence of the energy identity (|3.7ip , we can establish our fundamental energy inequality. 

Proposition 3.4 (The energy estimate). Suppose that the bootstrap assumptions (2.29a\) and (2.29b)) 
hold with e > and rj>0 sufficiently small. Letting K 



II go II 4 

II 9o II o ' 



sup S(s) + l f V(s)ds<£(Q) + CK 2 [ e" s : £(s)ds + 0(y/7) sup £(s) for te[0,T}. (3. 

0<s<i 2 J Q J q 0<s<t 



73) 



Proof. Throughout the proof, we will rely on the a priori bounds of Section [2 in particular, we will 
often make use of Lemmas 12.31 12.61 12.91 and 12.101 



Step 1. The estimate for Guopf in {3. 72\) We claim that 



\Guo P {\<CK / e^ s £{s)ds 



(3.74) 



Note that 



(-d N q t )R 2 j\B e h\ 



<C 



o Jr 



j-dwqt) 
-d N q 



(-d N q)\d 6 h\ 



<c 



dNqt 



-d N q 



£{s)ds. 



In order to bound the term 



On qt 
d N q 



we need a decay estimate for the numerator |9/v<7t|. The Sobolcv 
embedding theory would yield the bound \dNqt\L°=> ^ ||<Zt||2+<5 f° r 6>0, but by definition of our decay 
norm E/3, it is only the H 2 (tt) -norm of q t for which we have the desired decay. Thus, we arrive at 
the decay estimate for q t by using a comparison principle together with Theorem 1 in Oddson |40j ; 
indeed, in Appendix [B] we prove that 



\d N qt\ L °°<K 2 c x e-W 2 
It then follows from the bootstrap assumption (|2.29bl) that 

dNqt(s) 



(3.75) 



< ° A Cie , . < CK 2 e^ s , 



c\e 



-(A 1 +r,/2)s 



-dNq(s) 

which, in turn, establishes (|3.74|) . 

Step 2. Estimates for IZj , IZj, and Qj in \3. 71% Our objective will be to show that 



JQ 



o Jr 



<0(y/l) sup £(s) + S I Vis), forj = 0,...3. 

0<s<t 



(3.76) 



We establish (|3.76p for the most difficult case, j = 0. The case when j = l, 2, or 3 can then be 
proven in a similar fashion. The proof for j — is divided into three parts, and we shall begin with 
the term TZq. 
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Estimates for the integral J n TZq . As derived in (|A.9|) , the term IZq can be written as 



1=1 



fto := ME c ^d 6 - ^SV + ( W . fe A^), s ^ 

r-r * — , .. ' * .. — . ' <■ „ ' 

6 / 



= :/2 



= :Ii 



-(/i4), fe 9 6 g9 6 w i -/i^{9 6 ,a fe }g9 6 w i -/iV 

v y v J * * 



1=1 



cttfAfdP-yk (B 6 q + B 6 ^-v) 



+ d l d 6 ~ l wd l v (d 6 q + d 6 y-v)-d 6 ^-v t (d 6 q + d 6 ^-v) 



=:It> 



Estimate of JqIi- For the extremal case 1 = 5, 

B 5 A k Bq, k B 6 v l <\\d 5 Al\\ Li \\dq, k \\ L 4d & v% 

! 

<||f-Id|| 6 . 5 ||% jfe ||o. 5 ||9V|| 



< 



< l|g|U flZ2.p1/2 

<%e~^e£ + SV, 
o 

where we have used Holder's inequality and the Sobolev embedding theorem, as well as Young's 
inequality together with Lemma 12.31 for the last inequality. 

If 2 = 4, then Lemmas 12.61 and 12.101 and Corollary 12.111 show that 

a 4 A fe 5V.aV| < ||a 4 4llo||9 2 g lfe |Uoo||9V||o<|/ l | 4 .5||?IU. 5 © 1/2 

< e(\\q\\i 5 + V)<e{ee~'> t £ + V)<e 2 e~^£ + eV. 
The case when 1 = 1,2 or 3 are estimated in the same way and yield the same bound. 
Estimates of J a I k for A: = 2, 3, 4, 5. The following estimate holds: 

h + h + h + h ^te-^E + SV. 
i o 

For the integral of I2, an application of an L°°-L 2 -L 2 Holder's inequality together with Lemmas 12.31 
and 12.61 leads to 



< 



< 



( M , k AlA$), s \\ L ~\\d»* r \\ \\d b v% 



gl/2 

\fiAtA k ;\\ w i^\\q\\ 3 —^ 



£,1/2 < 



c 



-jt 



e£ + 5V. 



The estimates for terms I3, I4, 1$, and I§ are established in the same manner. Note that the 
commutator {B 6 ,A k }q, k in I3 is defined in p. 47b ) and has at most five derivatives acting on q, k ; 
moreover, the expression {B 6 ,d k }f = B 6 d k f — 8 k B e f is of the form J2i<\ a \<6 aa Baf ', where the a a 
are smooth uniformly bounded functions on the set a; = {x 6 f2| | < |ar| < 1}. 

Estimate Jn^7- We first consider the case that 1 = 6, and write 

r B e A'yv,l(B 6 q+3 6 ^-v)= [ B 6 A k v,lB 6 q + [ B 6 A k v,lB 6 ^-v. 
n Jn Jn 

- y ' v . ' 

Jl -h 
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Thanks to Lemma [3.31 we see that J\ < ||<9 5 A|jo.5||-Dv<9 6 <7||o.5- By linear interpolation and the 
Sobolev embedding theorem, ||-Df d 6 g||o.5 < IMUIklU + IMI2.5II9II6.5 ^ IMUIklle.s- It thus follows that 

■/i<|fe|6|| 5 . 5 ||^||3||g|l6. 5 < T NglHI§ + ^|g||g. 5 < 5 ^ {t) <SV + ee-^S + S(ee-^S+V), 

for some positive constant 7 > 0, where we have employed Lemmas 12.31 and 12.101 with Corollarv l2.11l 
As for the integral of J 2 , we again use Lemma 13.31 to deduce that 

where 7 > is given by Lemma 12.31 Now for the case that I = 5 in the integral of the term I7 , it 
follows that 



d a A?dv,i(d q+d°*-v) <\\d A«\\ L 4dv; k \\ L 4d b q +d^-v\\o 

cl/2 E l/2 e -0t/2 

< \\B 5 A%. 5 \\dv,l \\o. 5 \\d% + d e * -v\\o < —^Wvh.sE 1/2 < p t £ < y/ie-*£, 

where we used Lemma 12.31 again and the fact that (by definition off), \\d 6 q + d 6 ^ ■ v\\% < £. Hereby 
we used the estimate (|2.30p . The remaining cases 1 = 1,2,3,4 follow analogously and the estimates 
rely on a systematic use of Lemmas 12.31 12.61 12.101 and Corollary 12. Ill 

Estimate of J^-^s- For the case that I = 1 or 2, we have that 
■d l v(d 6 q + d 6 ^) <\\8 6 - l w\\o\\d l v\\ L ~,\\d 6 q 



d 6 - l wd l v(d 6 a + d e y) <\\d 6 - l w\\o\\d l v\\ L ^\\d e q + d 6 ^-v\\o< , ^^ El^e-^S 1 / 2 



<^^£+sv, 





while for the case that Z = 3, 4, 5 or 6, 

8 6 - l wd l v(d 6 q + d 6 ^-v) <\\d 6 - l w\\ L ^\\d l v\\ Q \\d 6 q + d 6 ^-v\\ <eV, 

! 

where we used the Sobolev embedding H 1+s L°° and Lemma |2"1)1 
Estimate of J^g- We see that 

8 e ^-v t (d 6 q +d 6 ^-v) <\\^\\4v t \\ L ^\\(^q+^-v)\\ <-^^El /2 e -' 3t / 2 f 1 / 2 <^ e - 7t/2 £:, 



<ee-~t t £ + 8V, 



with the decay rate 7 > given in Lemma 
Estimate of J^lZi. In the same manner, we find that JqT^i 



Estimate of the boundary integral J r Qq . We begin with the formula (|A.10[) (whereby we recall (|1.8|) 
n = AN=y/R? + R*n). 

5 



g = -d N qd 6 V ■nd 6 V-fi t -d N q 



dt 



Rd 6 h(- Rj + J2 c id a hd 6 - a ^ (h£- h g T)) 



a=0 



K, 



+ d N q—\{- Rj + ^c J a d a hd & - a (, - {hi~h e T)) 2l \ + J2ai(-d N q)d 6 ^-nd 6 - l (v-w)-d l n. 



1=1 



Estimate of J n K\ . Note that 
d N q t d 6 y-nd 6 if-h t 



K 4 



<\d N qt\L^\n t \ L ^\d 6 ^\l<\d N q t \ 1 \h t \ 2 ^<E -^£<ee-'' t £, 
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where we used the trace theorem and Lemma 



Estimates of J n K.% and J n A'3 . These two integrals are lower-order and thanks to Lemmas 12.6 
and 12.31 arc bounded by ee~ Jt £ + ST>. Note that |J| =l + 0(e) remains close to 1 due to the a priori 
smallncss bounds from Lemma 



Estimate of J n K^. The estimate of JqK^ requires some explanation, as it has the largest derivative 
count in Qq. In Appendix wc derive the identity 



d 6 ^ • n = Rjd 6 h- R,, + s ^c J a d a hd b ~ a £, ■ (h^-hgr), 



(3.77) 



where we recall that r is the unit tangent defined by (|1.3|) and Rj = RJ 1 . Substitution of (13.771) 
in the integral / r (— dNq)d 6 fy ■ hd 6 ~ l (v — w) ■ d l n then yields 



(-d N q)d b ^-hd b ' l (v-w)-d l h 



< 



(-d N q)8P- l (v-w)-&n 



(-d N q)0(d b h) ■ nd 6 ~ l {v -w)-d l n 



1~ 



{-d N q)Rd 6 hd 6 - l (v-w)-d l h . (3.78) 



The first and the second integrals on the right-hand side of (|3 . T8[) are easily estimated using Holder's 
inequality and the Sobolev embedding theorem, while the third integral on the right-hand side 
of (|3.78j) requires some care due to the presence of d 6 h. If I = 1 or I = 2, then 

(-d N q)RB 6 hd 6 - l (v-w)-d l h <|^9^5 6 / l | |\/^^ii|L-(|5 4 «| 1 + |5 4 / lt | 1 )|a z n| L o O 



<£ 1 / 2 || g ||J /2 (e e ^ t / 2 £ 1 / 2 +I? 1 / 2 )e< e 2 P(£,^) e -^ + <52?, 

where we have used Corollary 12. 11[ Lemma T2.61 and then Young's inequality for the last estimate. 
The case that 1 = 3, 4, or 5 follows similarly from Lemmas 12.31 12.61 and l2.10l The case 1 = 6 appears 
problematic because of the term d e n-r which, modulo coefficients, is essentially d 7 h, one derivative 
more than appears in £. The integral is, however, easily estimated thanks to the presence of an 
exact derivative, formed from the integrand d 7 h d 6 h. 

We set Jh = \J R 2 + and write the unit tangent to T(t) as t= J f ^ 1 (RT + hgN). A simple com- 
putation shows that 

n e = J- 2 (R 2 + 2hj + Rh eg )t. 

Since v — w = t- (v — w)i on T, we see that d 6 n ■ (v — w) = t- (v — w) B 6 n-i. We then write 

d 6 n- (v — w) = gid 7 h + 92, 

where g\ = t- (v — w) J^ 2 R, and where (72 is a lower-order term in v — w and has at most six tangential 
derivatives on h. We then write 

{-d N q)Rd 6 h(v-w)-5 6 n= [ (-d N q)R 9l d 6 hd 7 h+ [ {-d N q)Rd 6 hg 2 



T 



d[(-d N q)R gi ]\d b h\ 
Arguing in a similar fashion as for the case that 1 = 1 or 2, we see that 



(-d N q)Rd 6 hg 2 



{-d N q)Rd b h(v~w)-d b n 



<V~ee- 7t £. 



Step 3. Estimates for Sj, Sj, Hj in \3. 7i| ). We next prove that 



Jn 



Jr 



<0(Ve) sup S(s) + 6 I V(s), j = 1,2,3. (3.79) 

0<s<* 
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We will analyze the case that j = 1, as the estimates for the case that j = 2 or 3 follow in the same 
manner. We begin with the definition of Si given in (|A.23[) as 

5i:= c ahl id a d h t A 1 ld 5 - a dl-\, k d 5 v i + S l x 



0<o + i,<6 
a<5,6<l 



J2 dind 5 - l ^ t ■ 8 l v(8 5 q t + d 5 ^f v) - ^ c xi id l A\d h - l v\ h (8 5 q t + 8 5 ^ 



t-v 



(1X231) 



i=i 



i=i 



where S[ is a lower-order term given by 

S[ = ( nk AiA k r ), s 8 b %d 5 v l + {8 5 d t ,A k }q, k 8 5 v l + {d 5 Al l }q, k 8 5 v*+[xAjd 5 %A k r q, k {d\d s }v* 

-(flA^^d^qtd'v'+^d^diqid^d^ + ^id^dkjdtqd 5 ^ 

Most of the estimates arc completely standard and we focus on the more problematic terms, char- 
acterized by the highest number of derivatives applied to two out of the three terms in our cubic 
integrands. For illustration, in the first term on the right-hand side of (|A.23I) we analyze the cases 
(6 = 0,a = l) and (6 = 0, a = 5). If (6 = 0, a = 1) then we first integrate-by-parts and an L°°-L 2 -L 2 
Holder's inequality to find that 

r B 2 A k d 3 d t q^d 5 v l + [ 8A k 8 3 d t q k 8 6 v 



dA k d 4 d t q. k d b v 

Jn 

< \\8A k \\ wl ,~ \\d 3 d t q, k \\o(\\d 5 v l \\ + \\d 6 v\\o) 

< \h\ 3 .4 qt \\ 4 (\\d 5 v l \\ + \\8 6 v\\o)<eV, 

where Lemmas 12.31 and 12.101 have been used. If (b = 0,a = 5) then 



< \h\ 6 \\q t h\\8 5 v\\o < J^W^e-^^i/s < *±1 E + 6V , 



where we used Lemmas 12.31 12.61 and 12.91 The remaining estimates in the expressions (|A.23j) 
and (|A.26[) for S\ and S\ follow in the identical manner. As to the boundary integral of "Hi, 
we state the formula for the integrand derived in (|A.24[) as 



Hi := 28 N qd 5 h t Rj ^ d a h t d 5 ~ a £ -n + J^ a l (-8 N q)8 5 ^ ■ hd 5 v ■ fid 5 ~ l (v -w)-& 



(1X241) 



a=0 



1=1 



We consider the boundary integral of the first term on the right-hand side. We begin with the 
interpolation bound 

2,1/4 



^4<Nr Ns ^^)T/4-> 



(3.80) 



where we have used (]2.46[) to bound |9 3 /ij| and the definition of T> given in (|1.12[) . If a — A in the 
first term of the right-hand side of (jA.24[) . then 



d N q8 5 h t Rjd 4 h t dt-fi 



d N qd{\d*h t \')RM-h 



< 



d(d N qd£-n)\dV 



<l 



1\W 



pi/2 



< e 3/2 e - 7 *2,l/2 



< e 3/2 e -2 7t + £ 3/22, i 



where we have once again used Lemma 12.31 in second inequality, the estimate (|3.80[) . and Young's 
inequality. If ae{0,l,2,3} then 



8 N qd 5 h t Rjd a h t d 5 ~ a Z-n 



< \d N q\ L ~\8 5 h t \ \Rj\ L c°\8 a ht\o\8 5 - a Z-n\ L o 
2?i/2 



< 



19113 



x(t )l/2 



6 < e 3/2 e -^t v l/2 < e 3/2 e -2 7t + e 3/2 v 
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where we used Lemmas 12.61 and 12.31 and the same idea as above. The estimates for the second term 
on the right-hand side of (|A.24|) follow in an analogous vein, relying crucially on Lemmas 12 .61 and [2.31 
This finishes the proof of p. 791) . 

Step 4- The proof of the lemma is a direct consequence of the bounds (|3.74[) . (|3.76j) . and (|3.79|) . □ 

4. Existence for all time t > and nonlinear stability 

4.1. Structure of the proof. The basic goal in our strategy for global-in-time existence and decay 
of the temperature function is to prove that on any time- interval on which the bootstrap assumptions 
(|2.29a|) and (|2.29b|) arc valid, wc have that 

sup £{s)+ [ V{s)ds<C K £(0), 

0<s<t Jo 

where Ck > is some explicit constant depending on K, Upon choosing the initial data (qo,ho) suf- 
ficiently small, wc can obtain an improvement of the first bootstrap bound in (|2.29a[) . In Section |4~21 
we show the improvement of the bootstrap assumption on Ep in (|2.29a[) and in Corollary 12.141 wc 
have already shown the improvement of the bootstrap assumption (|2.29bj) . By a continuity argument 
this leads to a global existence result. 

In order to implement the above strategy, we start with the basic energy inequality given by (|3.73|) . 
Note however the presence of an exponentially growing term CK 2 J Q e vs £(s)ds on the right-hand side 
of (|3.73|) . That term appears by treating the terms / r (— dNqt)Rj\B 6 ~ 2j <9;? h\ 2 d6, j = 0,l,2,3 as error 
terms. By applying a straightforward Gronwall-type argument, this will be enough to guarantee 
that solutions to the classical Stefan problem p. 91) exist on a sufficiently long time- interval [0,Tr-], 
where the time Tk may be larger than the time of existence guaranteed by our local wcll-poscdncss 
theorem in |29j . As we explained in the introduction, by a sufficiently long time- interval, we mean 
a time Tk after which the dynamics of the Stefan problem (|1.9a[) are, in fact, dominated by the 
projection of the solution onto the first eigenfunction <pi of the Dirichlct-Laplacian. 

To prove global existence we need, however, more refined estimates that will show that the 
J r (—dNqt)R 2 j\d 6 ~ 2 -' h\ 2 d9 are in fact sign-definite for t>Tx, leading to the elimination of the 
exponentially-in-time growing bounds. First, in Section 14.31 we prove strict positivity of the term 
9jv<Zt at time Tk- Finally in Section |4.4[ we use a comparison principle to show that dNQt remains 
positive after time Tk- This allows us, in turn, to prove the uniform-in-time energy bound and 
extend the solution for all time t > 0. 

4.2. Boundedness of Ep. The following lemma shows that under the bootstrap assumptions, the 
bound on Ep + J^D^ds from (|2.29b[) is improved. 

Lemma 4.1. There exists a constant C and e>0 sufficiently small, such that if the bootstrap 
assumptions \2.29a}) and \2. 29b)) hold with such e and C , then 

Ep(t) + J D(s)ds<^Ep(0). 

Proof. We set 

x(t) = \\q(t)\\l + \\q t {t)\\ 2 + \\q u (t)\\ 2 and recall that D(t) = \\q{t)\\ 2 5 + \\q t {t)\\ 2 3 + \\q tt {t)\\ 2 1 . 
Step 1. Energy inequality for qu- From equation (|2.50b ). we see that 




where the forcing term fi is defined just below equation (|2.50|) . We next show that the right-hand 
side can be bounded by eD. We first focus on the term (^ t -v) u in the forcing function / 2 . Using 
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the product rule we obtain 

(*c»)iffe= / ^ttt-vqtt+ I 2^ t fV t q tt + I ^! t -v tt qu- 




For the integral A\ . we see that 



^ttt-vqtt 



<||*iit||o||w||L»||g«||o<|ftm|o.6|H|a||«ftt||o<ci5, 



where we used the bound (|2.54p to estimate |/ittt|o.5 by E 1 ^ 2 . The estimate \A2\ ^e£> follows analo- 
gously to the estimate for term Ai and the bound on ^3 follows from 



As 



<ll*t||i»l|utt||olktt||o<||*t||2<eA 



where we have used Lemma [2.61 to infer that || 1| 2 < e. All of the remaining terms in the forcing 
function fi can be estimated by a straightforward application of the Sobolev embedding theorem 
together with Lemma 12.61 (to guarantee the smallness of various Sobolev norms applied to the 
coefficient matrix (A^)k,i=i,2)- Thus, in summary, 

lj t hu\\ 2 + \\V*qu\\ 2 o<CeD. (4.81) 

Step 2. Elliptic estimates. We next prove that the quantities x and y are respectively controlled 
by II qtt || an d || V^^tt || q. Using the elliptic regularity estimate of Lemma 12^8| the elliptic equa- 
tions (|2.50p . and Lemma 121)1 it follows that 

IM|2<|M|o + ||/i||o, (4.82) 

and 

IM|4<||gJ 2 + ||/o||2<IMIo + ||/i||o + ||/o|j2. (4.83) 
A straightforward application of the Sobolev embedding theorem together with Lemma 12.61 implies 
that 

H/i|lo + ll/o|l!<ecW- (4-84) 

Hence, with (fQ2 )| -([£gg |) . 

x(t)<\\qtt\\o + ex(t), 
so that for e > taken sufficiently small, 

x(t)<htt(t)\\ 2 . 

Since 1 1 fx \ \ \ + 1 1 f \ \ \ < eD {t ) , the same argument provides 

D<\\qtt\\\<\\Vq tt \\l<\\Vmtt\\l, (4.85) 
the last inequality following from the uniform lower-bound of the matrix AA T . 

Step 3. Poincare inequality. The following bound holds: 

(Ax — 0(e))||y||g< ||V*/||g, (4.86) 

where V*=A T V and feH^fl). To see P~g6]) . note that the inequalities §ZM\) and (j!Q2"]) show 
that 

p-Id|| L ~<e, 

from which it follows that A^A?^ > (1 -0(e))|£| 2 for all £eR 2 . The Poincare inequality Ai||/||g< 
||V/||q for all qEH^tt) then concludes the proof. 

Step 4- The differential inequality and decay. From (|4.81[) and (|4.85l) we obtain that 

lj t \\qu\\l + (l-O(e))\\V^q tt \\l<0. 
Using the Poincare inequality (|4.86[) . it follows that 

j t \\qu\\ 2 + (2X 1 -O(e))\\qu\\l<0. 
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From this differential inequality, we immediately infer the bound 

lk tt Wllg<llto(o)|| e - (2Al -° (e)) *. 

From the elliptic estimate in Step 2, it finally follows that 

Since Ep(t) = x(i)e /3t and /3 = 2Ai — ?] < 2Ai — O(e) for e sufficiently small, it is now clear that we can 
choose C so that on the time interval of validity of bootstrap assumptions (|2.29a[) and (|2.29b[) we 
actually have the improved bound Ep{t) < ^e - ^*. □ 

4.3. Pointwise positivity of &Nqt at time Tk = C\iiK. 

Lemma 4.2. Assume that the solution (q,h) to the Stefan problem \1.9)) exists on a given time 
interval [0,T]. Let the bootstrap assumptions 112.29a]) and \2. 29b\) hold on that time interval with 
e > sufficiently small, and assume the smallness assumption \1.22\) for the initial data. There exists 
a universal constant C such that if T>Tk '■= ChxK , then 

-q t (T K ,x) > Ccie- AlT ^i(.T), xe Bi(0), 

where ipi is the first eigenfunction of the Dirichlet-Laplacian on Q and c\ = J^qoipidx. As a conse- 
quence, 

inf d N q t (T Kl x)>0. 

Proof. Step 1. Hardy-type estimate. As a consequence of the higher-order Hardy inequality (see 
Lemma 1 in [H]) and the Sobolev embedding theorem, for any /e7J 2 25 ( J B 1 (0))ni?o(Bi(0)), 

<C||/|| 2 . 25J (4.87) 



sup 

xG-Bi(O) 



where <pi be the first eigenfunction of the Dirichlet-Laplacian on the unit ball. 
Step 2. The Duhamel formula. Let 

oo 

i=i 

be the eigenvector decomposition of the initial datum qa with respect to the L 2 orthonormal ba- 
sis {ipi,(f2,...} associated with the Dirichlet-Laplacian on the unit disk -Bi(O). Writing the time- 
differentiated Stefan problem as a perturbation of the linear heat equation, we see that in fi, q t 
satisfies 

q tt -Aq t = N(q,h), (4.88) 

where 1 

N(q,h) := (ay -5ij)q t ,ij +h<lt,i + aij,tq,ij +b ltt qi+A%q, k w l +A%q, k w\ , (4.89) 
and the coefficients ay, bi are defined in (|2.70|) . Note that at time t = 0, q t (0) = Aq + \7q -wo; 
moreover, since Aipj = —Xj(fj and e tA is a linear semi-group , the Duhamel principle implies that 
the solution q t to (|4.88|) can be written as 

oo „ t 

-g t = c 1 Aie- Xl Vi + y]c 3 A J . e -^V 3 -e tA (Vgo-^o)- / e^ A N(q,h). 
T^i * ' Jo 

V ' 7~> 

We first prove that X(t) > for times t = C\nK, where C denotes a universal constant. We shall 
then show that at time t = C\nK, \Y(t)\ + \Z(t)\ is bounded by a small fraction of X(t). 

Step 3. Estimate of X. We begin by writing X as 

X(i, a; ) = c 1 A 1 e-^Vi(a:) + ciA 1 e- Al ViW(E^e^-^) t ^44). (4.90) 

V^— 'ciAi <i?i(x)/ 
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Our goal is to prove that the term 

a:= y^h e (x 1 -x j )t^M (491) 

'H'ciXi ipi (x) 

j — ^ 

is small. By Corollary 12.21 

IC'I 

^ < K for all integers j > 2 . (4.92) 

Cl 

Furthermore, using the normalization ||</jj||o = l, and the eigenvalue problem, Aifj=— Xjtpj, elliptic 
regularity shows that \\tp\\ 2 < Xj and that \\ip\\i < A 2 ; hence, linear interpolation provides us with the 
inequality 

H^ll2.25<A)- 25 . (4.93) 
Using (14.921) and (|4.931 together with the bound (|4.87p . we see that 

|£r|<CA-^A^ 25 e( Al -^)*. 

Since Ai < A2 < A3 < . . . , there exists a constant c* , uniform in j > 2, such that X\/Xj < (1 — 2c*). This 
implies that 

(Ai — Xj) < —2c*Xj for integers j> 2. 

In particular, for t> Chi-ftT 

A 2 25 1 

Cii:A 2 - 25 e- c * Ajt < CKXf 25 K- Cc * x > = C— < - 

for C chosen sufficiently large, but independent of K. (Recall that K > 1 since K > jj^jj* > 1 + Ai > 6.) 
Hence, from ()4.91j) and the previous inequality it follows 

-.00 -.00 -. 

i=2 j=2 

Plugging this into (|4.90|) . we obtain for any xG -Bi(O) 

X(t,x)>^c 1 X 1 e~ Xlt ip 1 {x)>0 7 t>ChiK. (4.94) 

Step 4- Estimates ofY and Z. The term Y satisfies the estimate 

l|y(i^)llL-<l|e tA (V go -wo)||2<e- Al 1Vg ^o||2e- Alt !| ( 7o!l3|^o!|2<ec 1 e- Alt , 

where we used the Sobolev embedding theorem together with the bound 1 1 go 1 1 3 ^ Kc\, which follows 
from 1 1 t/o 1 1 4 / 1 1 1 1 < K. Thus |Y(t,a;)| < \\X (t,x)\ with e sufficiently small. Next, to estimate Z 
which vanishes at the boundary, we have that 

Jo PiW Jo 



\N(q,h)(s)\\ 2 . 25 ds<Vt( / \\N(q,h)(s)g, 25 ds) 



1/2 



< 

In the above chain of inequalities, we have used the bound (|4.87|) for the second inequality, and the 
fact that \\e tA \\ H s^ H s < 1. 

We shall conclude our estimate by showing that 

f\\N{ q ,h){s)\\l 25 d S <c\ 1/b e^\ (4.95) 
Jo 

We recall that 

N{q,h) := (a lJ ~S. lJ )q t , l3 +b l q t ^ + a l: j tt q n: j +b lit q l + A%q, k W i + A\q, k w\ , (|4.89|) 
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and note that Z\ is the highest-order term with respect to the number of derivatives applied to q. 
Writing Z\ = (a — Id) D 2 q t , where Id denotes the identity matrix, we see that 

/ ll^l||i.25 ds ;SI|a- Id lli25lP 2 «t|l2.25^ SU P II a " Id II 2.25 / IMI4.25- 
JO 0<s<* Jo 

From the sharp estimate (|2.40p . we infer that sup 0<s<t \\a — Id||| 25~ c ii furthermore, for the term 

lktlU.25 we appl y the interpolation estimate \\qt\\i. 2 5 ^ ll&ll^lktll! 5 ^cj /5 e _ ^/ 10 ||g t ||^ 5 5 . 
Using Lemma ?2.W\ we then infer that 



JO 



the last inequality following from Holder's inequality and the fact that L \\qt\W 5 < e 2 by Lemma fe.lOl 
and the bootstrap assumption (|2.29a[) . 

Analogous estimates are applied to the term Z2 to finally deduce (|4.95[) . By (|4.95p and the above 
chain of estimates, it follows that 

Hence, at time T = C\nK 

1/10 

|Z(T,x)|<£V2 ln tfi/2 ClC Vio^ 

^ W^Tj20 ° 1/2]nKl/2 ^( x ) < \ci\e-^ T ^{x) < \x{t,x). 

Note that we have used the estimate c^ 10 < e^ 10 / F(K) 1 / 20 (which follows from ||g || <£{®) 1 ^ 2 and 
the smallness assumption (TTT^U)) as well as eo < e which is going to hold by our choice of eo- Observe 
that the very last inequality follows from (|4.94[) . The next-to-last bound is equivalent to 

£ < *! 

F(K)V W 4C 1 /2ln/v 1 /2xC*Ai' 

which then follows from the choice (|1.23|) of the function F(K) in Rcmark ll.3l The second inequality 
above follows from the estimate Ci<||g || < £(0) 1/2 < e/F(Kf/ 2 . 

Step 5. Finishing the proof. From the above estimates on X, Y , and Z it finally follows that for any 
ieBi(0), T = C\nK, 

-q t (T,x) > \X(T,x)\- \Y(T,x)\- \Z(t,x)\ > X(T,x) - ±X(T,x) - \x(T,x) = l -X{t,x) 
> c*c 1 e- XlT i Pl (x). 

Finally, since 9/v</Si>c for some c>0 uniformly over T and since tpi>0 in tt, it follows 
mi xer d N qt{T,x)>0. □ 

4.4. Proof of Theorem [TT2l Step 1. By Proposition ^. 41 with e sufficiently small, we conclude that 
sup £{t)+ [ V(t)dt<2£ {0) + CK 2 [ e nt £(t)dt, te[0,T], (4.96) 

0<t<T Jo JO 

where T is the maximal interval of existence on which the bootstrap assumptions ()2.29a|) and (|2.29bl) 
hold (with e sufficiently small). Our goal is to prove that on [0,T], the quantity £(t) is bounded 
from above by 2£(0)e 2CK *. We shall accomplish that by bounding £ (i) from above by the function 
g(t) :R + — >-K, which is defined as the solution of the differential equation 

g'{t) = CK 2 e^g{t), g(0) = 2£(0). 

Solving this differential equation, we obtain 

g{t) = 2£(0) e ^( e " t - 1 ) = 2£(0) e ( 1+ ^"" 1 * t / fc! ) Cif2t 
= 2£ (0)^ 1+o ^ CK2t < 2£(0)e 2CK2t , 
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where the convergence of the sum Y^ C k =2 11 l k 1 t k /k\ is guaranteed for times t<-^. Applying the 
integral Gronwall inequality to the difference £(t) — g(t), it follows from (|4.96[) and the previous 
inequality that 

£(t)<g(t)<2£(0)e 2CK2t 

for any t<T- Our goal is to prove that T>C\nK. Using (|4.96[) once again, we obtain the same 
smallness bound on J*D(s)ds to finally conclude that 



Jo 



sup £{s)+ [ V(s)ds<2£{0)e 2CK2t . (4.97) 

<s<t Jo 



0<s<t 

For t<C\nK, (|4.97|) and smallness assumption ()1.22|) on £(0) implies that 



sup £(s)+ [ V(s)ds<e/2. 

0<s<t Jo 



Moreover, by Lemma |4~T1 and since the bootstrap assumptions (|2.29a|) and (|2.29b|l are valid on [0,T] 
it follows that 



Ep{t)- 



J D(s)ds<^E (O). 



Thus, by the continuity of £ + Ep and the maximality of T, we conclude min{T,ClnK} = ChiK = 
Tk since the bootstrap assumptions are still satisfied at time t = ClnK (the argument is true as long 
as 77 above is chosen in such a way that -3= > C\nK). By the local well-posedness theorem from 

and the continuity of £ and Ep in time, we actually have the strict inequality T > Tk as we can 
extend the solution locally in time. We will argue by contradiction that T=oo. Assume T<oo. 

Step 2. Preserving the positivity of d^qt- We next show that d^qt >0 on the time interval [Tk,T[. 
This will be done with help of Lemma 14.21 and the maximum principle. We start by constructing a 
suitable comparison function, 

P(t,r) = Kl e~i Alt (<^(r) - k 2 (1 -r 2 )), (4.98) 
with positive constants ki,k 2 to be specified later. A straightforward calculation shows that 

3 1 

(a t -ay5„-& i a i )P = Kie"5 A i t [__A 1 ( / 3 1 -2K 2 Tr(a) 
3 

+ - Aik 2 (1 - r 2 ) - (a^ - 5ij )(fi-b- ( V(^i + 2k 2 x)] . (4.99) 

Observe that both (fi and (1 — r 2 ) vanish for r=l, the trace of the matrix a is very close to 2, i.e., 
an +a-22 = 2 + 0(e) and the coefficients hi are very small, i.e. |6| = 0(e). Note that the first and the 
second term in the parenthesis on the right-hand side of (|4.99[) are negative, while the fourth and 
the fifth term are small of order e. If r= \ x\ is close to 1, then the second term dominates the third 
term and if r is away from the boundary r = l, then one can choose /«2>0 so that the first term 
dominates the third term. It follows easily that there exists a k 2 > and some constant C\ such that 

{dt-aijdij-hdJP '<-CiKie-i Alt : . (4.100) 

It then follows from (|4.100p and ([2~7D1) that 

{dt-dijdij-bidiX-qt -P)> uhu +d t b l q l + d t A k l q lk w'+A^q, k wl) + C lKl e~^ t . (4.101) 

Note, however, that the term in parenthesis on the right-hand side above is a quadratic non-linearity 
and as such decays at least as fast as e _2/3 *: 

\\dtdij q,ij +d t h q l + d t A k ti q, k w l + A k q, k w\\\ L ™ 

< \\d t aijq,ij \\i +s + 110^^,111^ + \\d t A k l q, k w t \\ 1+ s + \\A k q, k w l t \\ 1+ s 

< C 2 E p (0) 1/2 ee~ 2(lt < C 2 c x ee- 2pt . (4.102) 

Now, using (|4.101[) and the above bound, we note that by choosing the constant Ki:=^-cie, we 
have that 

{d t - aijdij - bidi){-q t -P)> C2 Cl ee-i x ^- C 2 ciee- 2/3t > 0, 
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since 2/3 = 2Ai— rj> |Ai. The previous bound implies that —q t — P is a supersolution for the op- 
erator df — aijdij — bidi. Moreover, by the construction of P, we have —qt — P = on r — dBi(Q). 
Furthermore, at time Tr = Ch\K, we have by Lemma l4~2l and (|4.98[) . that 

(-q t -P)\ T= C lnK >Cc 1 e-^ T ^ 1 (x)~Cc 1 ee-^ lT Mx) + Cc 1 e K2 e-^ T (l-r 2 )>0 

for e sufficiently small. Thus, by Lemma \2. 121 there exists a constant m > such that 

-q t -P> m(l - r)e- (Al -° (e)) *, f > T K , 

or in other words 

-«t > m(l-r)e- (Al -°( e »* + Ccie(l-r)e-i Alt (^i^- K2 (l + r)) 

1 — r 

which readily gives the positivity of di^qt on the time- interval [Tk,T[ since the negatively-signed 
term in the parenthesis above has an exponential decaying weight in front of it. We conclude that the 
positivity of —q t at time Tk = ChxK is a property preserved by our bootstrap regime and moreover 
we get a quantitative lower bound on di^qt on the time interval [Tk,T[- 

Step 3. Conclusion. Thus for any t£ [Tr-,T[, the energy identity takes the form 
£(t) + f V(t) + \Y1 f I d N q t R 2 ,\d^dlh\ 2 dx 

J T K j = Q J Tk J V 

= £(T K ) + J2 I I {1U+Si}+^2 f [ {Ki+Si}+Yl f f {Gi+Ui}, 

where we formally define 54 = 54 = ^4 = 0. In particular, by the energy estimates stated in (|3.76j) 
and (|3.79p the right-hand side of the above identity is bounded by 

£(T K ) + 0{^1) sup £{s) + (0(e) + S) f V{s)ds. 

T K <s<t Jt k 

Note here the absence of the exponentially growing term in the above bound as opposed to their 
presence in Proposition I3.4I This is due to the fact that terms j^, j v dNqtR 2 j\d & ~ 2: ' d{h\ 2 dx, j = 
0,1,2,3 are positive and no longer treated as error terms. By absorbing the small multiples of 
sup Tjf<s<t £ (s) and J^, T>(s)ds into the left-hand side and using the positivity of dwqt from Step 2, 
we obtain 

sup £(s)+ / V(s)ds<2£(T K )<8£{0)e 2CK2TK (4.103) 

T K <s<t Jt k 

by (|4.97j) . Finally, we choose e in the statement of Theorem 11.21 so that e„<e 2 /2. Bound (|4.103|) 
and the condition £(0)<el/F(K) (with F{K) given as in ([T23")) ) imply 

f* e 2 
sup £{s) + / V(s)ds<—. 

T K <s<t Jt k 2 

Together with Lemma 14.11 and Corollary 12.141 we infer that the bootstrap assumptions (|2.29al) 
and ()2.29b[) are improved. Since £ is continuous in time, we can extend the solution by the local 
wcll-poscdncss theory to an interval [0,T +T*] for some small positive time T*. This however 
contradicts the maximality of T and hence T=oo. This concludes the proof of the main theorem. 

5. The ^-dimensional case on general near-spherical domains 

In this section we briefly sketch the set-up of the problem in general dimensions and explain 
how to adapt the arguments from the 2-D case to the 3-D case. Let f2(t)cR d be an open simply 
connected subset of R d , d> 2. The moving boundary T(t) = <9fi(i) is parametrized as a graph over 
the unit sphere § d ~ J 

T(t) = {x\ x = R(t,£)£=(l + h(t,t))t, eeS"- 1 }. 
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Initially Rq(£) is assumed to be close to 1, i.e. Ro(0 — 1 = ho{0 = 0(e). We shall assume that 
Do is diffeomorphic to i?i(0), where $:f2 — ^ f^o is the diffeomorphism mapping of the unit ball 
onto the initial domain. Moreover, let W denote the family of diffcomorphisms from the initial 
domain flo to the moving domain fl(t), satisfying the harmonic equation A\& = and the boundary 
condition ^(To) = T(i). We shall pull back the Stefan problem onto the unit ball £?i(0) via the map 
: Bi(0) — > Q(t) given as a composition of ^ and $: 

Upon defining q, v, w, and A just as in Section (|1.9p . the Stefan problem (|1.1[) takes exactly the 
same form as (|1.9p . Abusing the notation, the normal velocity V(T(t)) is now given by 

V(T(t)) = 



Here |Vi?|g d _i stands for the squared norm of the Riemannian gradient of R(t,-) on the unit sphere 
§ d_1 , which is a coordinate invariant expression. The gauge equation for ^> transforms into 

A 4 -i*=o, v(t,§ d - 1 ) = r(t) 

due to the assumption A^ = and the definition of fy. This easily implies the optimal trace 
bound H^Hijam) < |$|j/»-o.B(f2) due to the the smoothness of $ and the closeness assumption 
— Id||/fs <e, with s sufficiently large. When d — 3, the Sobolcv embedding theorem requires 
us to raise the degree of spatial regularity in the definition our energy spaces by one derivative. 

The second key observation is that the lower bound for the quantity x{t) is obtained in the same 
way as in the case that d = 2, from Lemma \2. 121 We v\ denote the first eigenvalue of the operator 
— A$-i, which is the pull-back of the negative Laplacian from the initial domain Oo to Bi(0). By 
Lemma 12.121 we obtain that 

X(t)>c ie ~ xt , 

where |A— v\ \ < 0(\h — ho\ w 2 ^ + |/itU°° ) = 0(e). Since ||D4> — Id||^fs <e for s sufficiently large, we 
have \v\ — \\\ <e, where we recall that Ai stands for the first eigenvalue of the Dirichlet-Laplacian. 
Together, the two previous estimates imply the analogous conclusion of Corollary 12 .141 namely 

x{t)> Cl e-^- x ^\ Ai=0(e). 

Let 3 l denote the tangential component of B % restricted to S 2 . To each multi-index a= {0.1,0.2,0.3) 
we associate the tangential operator B a = B ai B a2 B a!i . With d = 3, we define 

SaD(t) = e aD (s,h)(t):= 

\ y, y 1/2 d s d b t v\\k+\ E \(-d Nq ) i/2 RjB s d b t h\i l + 1 - e w^id^u+B^-v)^ 

|q|+26<6 |5|+26<7 |5|+26<7 

E H(i-M) 1/3 ^aMii 2 +^ E m-^^q+dsa^-v)^ 

|<3|+2b<6 |a|+2h<7 

and 

-D3D(t)=V 3D (q,h)(t):= 

E y 1/2 B s 8'v\\li+ E \(-B N q) 1/2 R,jB s B b t h t \l l+ E y 1/2 (B s B b qt + B s 8^fv)\\ 2 Ll 

|5|+26<7 |S|+2b<6 |5|+2h<6 

+ E \\(l-ri 1/2 d S B b t v\\l l + E Ul-^iBaB'qt + BsB^fv)]^. 

|S|+2b<7 |<S|+2b<6 

The lemmas of Section [2] carry through analogously, as do the energy estimates of Section [3J By the 
continuity argument of Section 01 we arrive at the 3-D version of our main theorem: 

Theorem 5.1 (The 3-D case). Let {qo,ho) satisfy the Taylor sign condition U.18\) , the strict positiv- 
ity assumption \l.l r ty , and the corresponding compatibility conditions. Let 1 1 (?o 1 1 4 / 1 1 <?o 1 1 < K. Then 
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there exists an eo = £o(K) > and 5q > such that if £(qo,ho) < eg, then there exists a unique global 
solution to problem il.9\) , satisfying 

sup s 3D (q(t)Mt))<csl 

0<t<oo 

for some universal constant C>0, and 1 1 Q 1 1 ^5 (^b x (o) ) — C e_,3 *i where /3 = 2Ai — Ceo and Ai is the 
smallest eigenvalue of the Dirichlet-Laplacian on the unit ball B(0,l)cl 3 . The moving boundary 
F(t) settles asymptotically to some nearby steady surface V and we have uniform-in-time estimate 

sup \h-h \ 5 . 5 < V^o 

0<t<oo 

Appendix A. Proof of Proposition 13. II 



To prove the energy identity of Proposition 13. R we start by applying the differential operator 
of the form d^^dl to the equation (jl.9bl) . For j = 0,l,2,3 we multiply it then by 8 6 ~ 2: >dl and 
integrate-by-parts. Additionally, if j = 1,2,3 we apply the operator 8 7 ~ 2: >d J t to (|1.9bj) . multiply by 
d 1 ^ 2j d\~ v l , and again integrate-by-parts. 

Based on these two cases we distinguish between the two different types of identities. 

A.l. Identities of the first type. Recall that /i : f2 — » R is a C°° cut-off function with the property 

fx(x) = if \x\<l/2; (jl(x) = 1 if 3/4 < \x\ < 1. 

Applying the tangential differential operator /id 6 to the equation (|1.9b|) . multiplying it by d 6 v z and 
integrating over H, we obtain 

5 

(/i«5V + ^A k q, k + ^A k 8\, k , V) L2 =J2°l {0A k d^ l q, k , flV) L2 , 

1=1 

where c; = ( j ) . Recalling (|2.47[) , we write 

d 6 A k = -A s l d 6 ^ s A k + {8 6 ,A k }, 
where {d 6 , A k } denotes the lower-order commutator defined in (|2.47b ). With this identity, we obtain 
( M d 6 A,^,aV) i2(n) = -^Atd & ^ s AU. k ,d & v l ) L2{n) + ^{d 6 ^}q, k ,d 6 v l ) L2{n) 

= - f q, k Atd 6 ^ r A k ;B 6 v 1 N s + f f iAtd e ^ r A k r q, k B e v l s + f Ti 
Jr Jn Jn 

= - f q^d 6 ^ r A k d 6 v l N s ^- [ f iA*d 6 y r v r 8 6 vl s + [ 71, (A.l) 
Jr Jn ' Jn 

where we have integratcd-by-parts with respect to X s for the second equality, and have used the 

identity v r = — A k q. k for the last equality; the error term 7i is given by 

71 = ( m , k AlA k r ), s 8 6 ^ r J e v l + fi{B e , A\ }q, k , <5V + [iAt8 6 * r A k ;q,k{d 6 ,d s }v\ 

Furthermore, integration-by-parts with respect to x k yields 



(flA k 8 6 q, k ,8 6 v l ) L2 = / iiA k d k d\d G v l + / [iA k {d e ,d k }q3 6 v l 

, _ M J n (A.2) 

= - / vA k 8 6 q8 6 v\ - / (M"U<9VV+ / nA k {d 6 ,d k }q8 6 v\ 



n 



where we have used d 6 q = on T, and where {d 6 ,d k } denotes the lower-order commutator. Sum- 
ming (jA.ip and (|A.2|) . we find that 



(fid 6 A k q tk + iiA k dS^d & v l ) L2[n) = -j q^ k Ai8^ r A k d & v l N s 

[iA k d C y k (d 6 q + d 6 y-v)+ f (Ti-^A^^qB^+fiA^d^dkjqd^v 1 ). 



The first two terms on the right-hand side of (|A.3[) will be the source of positive definite quadratic 
contributions to the energy. To extract the quadratic coercive contribution from the first integral 
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on the right-hand side of (|A.3|) . note that q,k = N k djs[q on T, and also recall from (|1.8[) the normal 
vector n = A T N. Thus 

-J q. k Ald 6 WA k ;B 6 v l N s = J {-d N q)5 6 y r h r d e v l h l = J (-3 N q)B 6 ^-hB 6 vh. (A.4) 

Using the boundary condition (jl.lOp , we reexpress B 6 v-n as 
B 6 v-h = B 6 w-h + 8 6 (v~w) - h 

6 

= 5 6 w-n+d 6 ((v-w)-n)-^2aid 6 - l (v-w)-d l h 

1=1 



=o 

6 



= dPVt-n-^2aid 6 - l (v-w)-d l 
i=i 

Due to the above identity and (|A.4j) . we obtain that 

f (-d N q)d 6 ^-fid 6 vfi= f (-d N q)d 6 ^-fi d 6 ^ t -n- Va, f {-d N q)d 6 ^-nd & -\v-w)-d l 
Jr Jr J=1 Jr 

= I (-d N q)~\BH-n\ 2 dx' + f dNqd^-nd^-ht-y^ai f (-d N q)B 6 V-hB 6 -\v-w)-B l 
Jr 2 at Jr ^—^ Jp 

(A.5) 

Recall that h = J- 1 (RN-R 8 T) = J-~ L (N + hN-h 8 T). Thus, using #(t,g)=N+h(t,£)N, we obtain 
via the Leibniz rule 

8 6 V-n = [B 6 N + B 6 (hN)]-[(l + h)N-h e T] J" 1 

5 

= ( -R + RB 6 h + c a B a hB & - a N ■ (hN - her) ) J" 1 

a=0 

5 

= {-R., + RjB 6 h + Y,cid a hB 6 ~ a N-(hN-h e T)), 

where we have used the relations B 2 N = —N and N ■ t = and also denoted c J a =c a J~ 1 (recall Rj = 
RJ -1 ). From here we obtain 

, 5 5 

+—[(-Rj + 2RjB 6 h+J2 c iB a hB 6 - a N-(hN- h e r)) (-Rj + ^ct B a hB 6 ~ a N ■ (hN-her))]. 

a=a a=Q 

(A.6) 

Thus, going back to (|A.5|) . we obtain 

J(-d N q)& t *-n8 B * t -n= ~J^-d N q)R 2 \3 6 h\ 2 + ^ ^8 N q t R 2 j\3 6 h\ 2 + ^T 2 , 
where the error term 7~2 is given by 

7 5 5 

T 2 = (-3 N q) — [ ( - Rj + 2R jB 6 h + ct B a hB 6 - a N ■ (hN - her)) ( - Rj + ^ c J a B a hB & ~ a N -{hN -her))}. 

a=0 a=0 

As to the second term on the right-hand side of (|A.3[) , note that 

6 6 

A$&v% = B 6 (A k v: k )-Y,ciB l A k B e -% = -B 6 (qt + v-w)~Y, c i BlA i B 



1=1 1=1 
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where A*v i k = — div*u = — (qt + v-w) by the parabolic equation (|1.9a[) . Thus 

- f h iA k B 6 vi k (B 6 q + B 6 ^-v) = [ f j,d 6 (q t + ^fv)(d 6 q + dH-v)+y2ci [ nB l A k B 6 ' l v l k (B 6 q + B 6 ^-v) 
Jn Jn i=1 Jn 

■ [ fi(B e q + B 6 ^-v) 2 + [ ^(y j d l d 6 - l -$ t -d l v-d 6 -i>-v t )(d 6 q + d 6 y-v) 
Jn Jn J=1 

Vc,/ fiB l A k B 6 - l v\(B 6 q + B 6 ^-v) 

, ; Jn 



ld_ 

2dt Jn 

6 



Combining (|A~3]| - (|A"7F|) we obtain 



n 2 at J r 2 at Jq 

~ [ (-d N q)\J2c a d a Rd 6 ~ a N-fi\ 2 dx' + I K Q + f Go 
Jr a=0 Jn Jr 



with the error terms TZq and Qq given by 
5 



(A.7) 



(A.8) 



K =/i^ Q 5'Afa 6 -' g , fc a 6 V I + { m , k AlA k r ), s d^ r B 6 v i + i i{B & ,A k i }q >k d e 
i=i 

+ h tA!B 6 ^ r A k ; q , k {B 6 ,d s }v l -( h iA k ), k B 6 q B 6 v l 

6 

- nA k {<9 6 , d k }qB 6 v l i c iB l A k B 6 - l vi k + did 6 ~ l w ■d l v-BH-v t ) (B 6 q + d e ^-v); (A.9) 

i=i 



Go = -d N q8 6 y -nd 6 y-n t + (-B N q) — \lljd 6 h( -R.J + J2 ciB a hB e ~ a N -(hN- hgr)) 



-dnq 



dt 



-Rj + J2 c t d a hB 6 ~ a N ■ (hN - hgr)) \+J2 ai (-dNq)B 6 ^ ■ fid 6 - 1 (v-w)-B 



a=0 



a=0 
6 



(A.10) 



1=1 



Let now a = (0:1,0:2) be an arbitrary multi-index of order 6. Applying the operator (1 — /j,)d a to (ll.9bl) 
and multiplying by d a v l , we obtain 

((1 - M )c>V + (1 - »)d a A k q, k +(1 - fi)A k d a q, k , d a v l ) L2{Q) 
= - J2 cp{(l-n)df'A$ff^f , q, k ,ePv i ) La 

0<fi<ct 

In the same way as above we arrive at the following energy identity 

1 d 



(1-H)\d a v\ 2 dx+-— / {l-Li){d a q t + 8 a ^-vydx= / TZ dx, 
2 dt j n 



(All) 



where 



Ro = (l-li) J2 c 8PA k d a -Pq, k B a v l + (1 - m) ((g,^^) !S S^dV + {3 a ,A-}q,kB a v 

0</3<a 

-(1-M) {cpd p J%d a - p v\ k + dpd a - p w-dPv-d a V-^){d a q+dP , *-v), 

0</3<a 



(A.12) 
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Summing the identities (|A.8|) and (|A.lip , with j = we arrive at 

f ^~ 23 d 3 v\ 2 + J2 I 0--^W^\ 2 + ~ I (-d N q)R 2 j\d e ^dih\ 2 

|a|=6-2j 

+ ~T I ^(p Q -^dlq + B 6 - 2l dl^-v) 2 + -^r- V / (l-^)(d a d 3 q + d a d 3 ^-v) 2 

(-d N q t )R 2 \d e - 23 d 3 h\ 2 + f (7^+7^-) + f G J7 (A.13) 

Jq Jr 

By imitating the same calculation as above we obtain the remaining error terms. With j = 1,2,3 the 
formulas for IZj, IZj, and Qj in (|A.13[) read 

Uj = J2 vd ab d a d b t A><;d 6 - 2: >- a dr b q,kd e - 2: >d 3 t v l + { nk A°A k ), s d 6 - 23 d 3 y r 8 e - 23 d 3 v l 

0<a+b<6-j 

+ nAid Q - 23 d 3 t V A% k {d 6 ~ 23 , d s }d 3 v l + ^{d 6 - 23 d 3 , A k }q, k d 6 - 23 d 3 v l 
- (Mf ),k d e ~ 23 d 3 qd 6 - 23 d 3 v l - [iA\{W- 23 ,d k }d 3 qd 6 - 23 d 3 v l 

-11 ]T (d ab d Q - 23 - a d 3 - 1 'A k d a d h t v% + d ab B a d^ t ■ B 6 - 23 ~ a dr b v - B 6 - 23 d 3 ^ ■ v t ) 

0<a+b<6-j 

x (d 6 - 23 d 3 q + d e - 2j d 3 t y-v); (A.14) 

5-2j 

Qj = -dNqd 6 - 23 ^^ ■ nd^ 2j d 3 t ^ ■ fit + d N qd t \d 6 ~ 23 d 3 hRj(-Rj + ^ d J a d a d 3 t hd & - 23 - a N ■ n) 

A 5 ~ 2j 



a=0 



+ rfah(-aAr g )a 6 - 2j a t J *-na 6 - 2j - Q 1 9r''(«-w)-5 a l 9 t f 'n. (A.15) 

0<a+fa<6-j 

0</J<a 

-(i-m) (c(i^r ,i i); l +(i-/i)(ip9"A-^i>4 <, ri) t )(a°9+3 a $-i)). (A.i6) 

0</3<a 

A. 2. Identities of the second type. Applying d 5 dt to (|1.9b|) and computing the L 2 (f2)-product 
with ^<9 5 v* we obtain 

{^vl+^Al^+nA^B^q^t^ 5 ^)^ = J2 Cab(vB a 8 h t AlB 5 - a a]- h q, k ,B 5 v l ) L2 , 



0<a+b<6 
a<5, fa<l 



where c ab are constants appearing due to the usage of Leibniz product rule above. Recalling (|2.47p . 
we write 

where {d 5 dt,A k } stands for the lower order commutator defined in (|2.47[) . With this identity, we 
obtain 

(/xd 5 At^,av) L2(nr -(M^ 5 ^ 

= - j q. k Aid b %A k r 8 5 v l N s + J^Atd b %A% k d b v] s + Jllx (A1?) 

= - / q, k AtB b %A k r d b v l N s - j nA s t d 5 V r t v r d 5 vi s + f U x , 
Jr Jn ' Jn 

where wc have integratcd-by-parts with respect to x s in the second equality and we have also used 

the identity v r = ~A k q, k to write the last line more concisely. The error term U\ is given by 
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Furthermore, integrating by parts with respect to x 



r T T r (A- 18) 

= - / »A*8 5 q t 8 5 V* k - / (^A k ), k B 5 q t B 5 v l + / ^A k B b q t {B b ,3 k }v l + / fiA k {B 5 ,d k }q t B 5 v\ 
Jn Jn Jn Jn 

where we have used B 5 q t =0 on T. Summing (|A.17[) and (|A.L8|) . we obtain 

(^A^^+M^^aV) L3(n) = - J^q. k Aid 5 y r t A k r 3\m s - J^A*d 5 v%(8 5 q t +d 5 * t -v) 
+ [ (wi -(Mf),^ 5 g t a^ ^ +Mfa 5 (Z^ {5^a fc }^; ^ +M ^ fe {a^9 fc }q ^ a 5 u ^ ). 

Jn 



(A.19) 

The first two terms on the right-hand side of (|A.19[) will be the source of positive definite quadratic 
contributions to the energy. To extract the quadratic coercive contribution from the first integral 
on the right-hand side of (|A.19|) . note that q, k = N k dwq on T. Thus 

-J q^AlB^A^^N 8 = J (-d N q)B 5 %fi r B 5 v l n l = J (-3 N q)B 5 ^ t -nB 5 vn. 

Just like in Section |A~T1 - as in the identities leading up to (|A.6|) - we obtain 

4 4 

{-3 N q)8 5 V t ■hd 5 vh=\d 5 h t \ 2 R 2 j + \Y^ c ,] a B a h t B b ~ a N -n\ 2 + 2B 5 h t Rj ^ c J a B a h t B^ a N ■ n 

a=0 a=0 

4 

+ a i {-B N q)B 5 ^ ■ fiB 5 v ■ fiB 5 - 1 (v-w)- B l h, (A.20) 



1=1 



where c J a = c a J 1 and c a are some universal constants. As to the second term on the right-hand side 
of (|A.19|) , note that 



A k B 5 v% = B"{A k v] k ) - ciB l A k B 5 - l v> k = -B 5 (q t +vw)-J2 ciB l A k B 5 - l v[ 
i=i i=i 

where A k v l k = — div^v = — (q t +v ■ w) by the parabolic equation (|1.9a[) . Thus 

- / LiA k B 5 v\(B 5 q t +B 5 *fv) 
Jn 

= / tlB 5 (q t + ^fv)(B 5 q t +B 5 ^fv)+y2ci [ vB l A k B 5 ~ l v] k {B 5 q t +B 5 ^ t -v) 
Jn J=1 Jn 

= f ^qt+B^fvj' + y^di [ fiB 5 - l ^fB l v(B 5 q t + B 5 ^fv) 
Jn ,_ n Jn 

5 



(A.21) 



+ Y,ci f fiB l A k B 5 - l v l k (B 5 q t + B 5 y t -v) 

1=1 Jn 



Combining (|A~19l) - (|A~2~T1) we obtain 



\— [ fx\B 5 v\ 2 dx+ [ (-B N q)\B 5 ^fh\ 2 dx , + [ tJ,{B 5 q t + B 5 ^fv) 2 dx 
2 at J n J r J n 

4 r r (A.22) 



f d N q\Yc J a B a h t B b - a t-n\ 2 + ( s x + \ui 

Jr a=0 Jn Jr 
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with the error terms Sj and ~Hj given by: 

S o = d ab iid a d b t Afd 7 - 2: >- a dr b q,kd 7 - 2: >dr 1 v l + (fiq, k A?A*), s B 7 ' 23 Bl^ r B 7 ~ 23 B^v 1 



0<a+b<6-j 
a<7-2j, b<j 



+ M {a 7 - 2 ^> J fe } g , fe ,d 7 - 2j a^ 
-{^), k d 7 - 2l diqB 7 ^dr 1 v l +^A'id 7 - 23 diq{d 7 - 2 ^d k }dr 1 v l 

+ Mf {B 7 ~ 2 \ B k }df qd 7 - 21 dr 1 v* - (<9 7 ~ 2j B\ q + d 7 - 2] d}V-v) 

x d ab (fid a d^t-d 7 ~ 2j ~ a dt 1 ^\ + fid 7 - 2j - a dt 1 ^ b Afd a d b v i k ). (A.23) 

0<a+6<6-j 

n 3 =2d N qd 7 ^dlhR,j c J ab B a d b t hd 7 - 2: >- a dl- b N-n 

0<a+b<7~j 

+ {-d N q)B 7 - 2l dl^-n J2 d ab B 7 - 2j - a d{- 1 - b (v-w)-d a a b n. (A.24) 

0<a+6<6-j 

Note that the first line of (|A~24)l appears as an expanded difference between two positive definite 
expressions (— 3Nq)\B 7 ~ 2j 3 J t * -«| 2 and (— dNq)\d 7 ~ 2j 8} h\ 2 . We do this just like after (|A.5[) using 
the formula n= J _1 (iV + hN — hgr) and the parametrization *&(t,£) = (l + h(t,£))N. Let now ct = 
(ot\,a2) be an arbitrary multi-index of order 5. Applying the operator (l — /j,)B a 3t to (|1.9b[) and 
multiplying by 3 a v l , we obtain 

((1 - H)d a v\ + (l-fi)d a Al t q, k +(l-^d a q, kt , d a v l ) L2{n) 

= - Y c^ b {{i- l j)d^d h t A k i d a -^dl- b q, k ,d a v i 



1 L- 



O<|0| + I><5 
P<a;b<l 



In the same way as above we arrive at the following energy identity 

\jt I ( 1 -V)\d a v\ 2 dx+ [ (l-fi)(d a q t + d a ^fv) 2 dx^ I S lC lx, (A.25) 

where 

Si= J2 C P»0-- »)dPd b t A%d a -Pdt- b q, k d a v i + (1- ^)({q^AiA k r ), s d a %d a v l + {d a d tl A k ;}q, k d a v' 



0<|3|+ b<5 
P<a;b<l 



(A.26) 



- (i - M ) Y c ^3 a ~^ t ■ B l v(B a q t +d a *fv)-(i-fi)Y cmd e A^d a -^v% (B a q t +B a ^ t - v) 

0</3<a 0<a 

For a general j G {1,2,3} we have 

Sj= Y d pb (l~^8 b A^8 a -P8r b q.kB a v i + (l^ f i)U 1 

O<|0| + b<7-2j 
0<a;b<j 

J] d pb fid fS d b ^fB a - 3{~ 1 ~ b v(B a 8{q + B a Bi^-v) ( A27 ) 

0<|/3|+fc<|a| 

-(1-A*) E dp^B^Br^A^Btv^Blq + B^l^-v). 

Q<\0\+b<\a\ 

Summing the identities (|A.22[) and (|A.25[) we arrive at 



[ v\3 7 - 2 >8lv\ 2 + Y [(l-t*Wdiv\ 2 }+ [(-d N q)R 2 j\d 7 -^dih\ 2 
+ f n(d 7 ~ 2j diq+d 7 - 2j di^-'"f+ Y ! (l-v)(d a 3 3 t q + B a 8 3 t y-v) 2 (A.28) 

|a|=7-2j 

= / (5,+^) + / Wj, i = 1,2,3. 
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Summing the identities (|A.13[) for j = 0,1,2,3 and (|A.28[) for j = 1,2,3, we conclude the proof of 
Proposition 13.11 

Appendix B. Proof of the inequality ([3.75 J) 

We use the comparison function P defined in (|4.98p with the same «2 and K\. Note that n\ = C*eci 
is defined as a multiple of c\ for some constant C* > 0. Using (|4.100p and upon possibly enlarging 
C*, we infer (d t — Oy — &»)(— qt +P) < 0. Theorem 1 from [40] guarantees 

-q t + P<C C 1 pe^ +C ^ t , (B.l) 

where p(r) = l — r stands for the distance function to the boundary T. Note that the constant «i 
in the definition (|4.98|) is chosen right after (|4.102[) . It is in particular proportional to Ep(Q) 1 / 2 < 
1 1 (7o 1 1 4 - By definition of K we have that ||golU < -^Ikollo- Since however ||9o||o < Kc\, we obtain 
P/p< CK 2 cie _3Al '/ 2 . Similarly, the constant Cq is proportional to the L°°-norm of the initial datum 
for —q t + P, wherefrom we again obtain Cq<K 2 c\ by the same argument as above. Dividing by p 
in (|B.1[) . from the above inequality, we infer that 

\d N q t \oo<CK 2 Cl e^ +c ^. 

This proves the inequality (|3.75l) . 
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